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Abstract 

We set up the singular initial value problem for quasilinear hyperbolic Fuch- 
sian systems of first order and establish an existence and uniqueness theory for this 
problem with smooth data and smooth coefficients (and with even lower regular- 
ity). We apply this theory in order to show the existence of smooth (generally not 
analytic) T 2 -symmetric solutions to the vacuum Einstein equations, which exhibit 
AVTD (asymptotically velocity term dominated) behavior in the neighborhood of 
their singularities and are polarized or half-polarized. 
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1 Introduction 

Fuchsian formulations have proven to be very useful for studying the behavior of cos- 
mological spacetimes in the neighborhood of their singularities. Introduced into general 
relativity almost fifteen years ago by Kichenassamy and Rendall [30] , these formulations 
have been used primarily as a tool for showing that within certain families of solutions 
of the Einstein equations (defined primarily by the invariance of each member of the 
family under a fixed isometry group), there is a large collection of solutions which ex- 
hibit AVTD (asymptotically velocity term dominated) behavior. Roughly speaking, a 
spacetime shows AVTD behavior if, in a neighborhood of its cosmological singularity, 
the evolution of the spacetime metric field of the solution approaches the evolution of 
a model metric field which (relative to some choice of spacetime coordinates) satisfies a 
system of ordinary differential equations (ODEs) deduced from the Einstein equations 
by suppressing spatial derivatives relative to time derivatives. The detection of AVTD 
behavior has proven to be a very useful step towards verifying that the strong cosmic 
censorship conjecture holds for certain families of solutions of the Einstein equations 
[18, 35]. 

Fuchsian formulations are effective in studying the possible presence of AVTD behav- 
ior since they are designed specifically to handle singular systems of partial differential 
equations (PDEs) or, equivalently, PDE systems in the neighborhood of their singular- 
ities. While it is not easy to identify the location of singularities in generic spacetime 
solutions of Einstein's equations, for certain isometry-defined families of solutions — e.g. 
spatially homogeneous solutions, Gowdy solutions, T 2 -symmetric solutions, and many 
families of [/(l)-symmetric solutions — one can use areal coordinates or special forms 
of harmonic coordinates to locate the singularities. If the Einstein equations are then 
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reduced relative to these symmetries and expressed in terms of these coordinates, then 
the resulting PDE system takes a singular form in the neighborhood of the singularity 
which is amenable to a Fuchsian formulation and analysis in the form of a singular initial 
value problem — presuming that various further conditions are met. 

Most of the earlier applications of Fuchsian formulations to families of solutions of 
Einstein's equations have presumed that the spacetimes are analytic [30, 24, 26, 4]. This 
is not surprising, since Fuchsian formulations for generic systems of PDEs were initially 
developed with analytic PDE systems in mind [29, 31]. It is important, however, to 
extend studies of AVTD behavior and strong cosmic censorship beyond analytic space- 
times and to consider if they also hold for spacetime solutions which are smooth, but 
not necessarily analytic. 

There are two sets of results (prior to our work) known to the authors concerning the 
existence of solutions to quasilinear Fuchsian PDEs in smooth or finite differentiability 
regularity classes. 1 In the first of these, proposed by Claudel and Newman [19], the main 
result is that if a number of quite restrictive technical conditions are satisfied by the 
PDE system, then the Cauchy problem is well-posed for data specified at the singular 
time. As noted in [34], these restrictive conditions are not generally satisfied by the 
PDE systems corresponding to the Einstein equations for the Gowdy spacetimes, for the 
T 2 -symmetric spacetimes, or for most other families of spacetimes under consideration; 
hence the Claudel and Newman results are not useful for our present purposes. 

The second set of results concerning smooth solutions of Fuchsian systems are those 
proven by Rendall. In [34], he develops a Fuchsian-based approach that is applicable 
to both semilinear and quasilinear equations, and he uses it to establish the existence 
of a class of smooth T 3 Gowdy spacetimes which exhibit AVTD behavior. In Rendall's 
approach, one performs a series of reduction steps in order to obtain a symmetric hy- 
perbolic system, and one then proves the existence of smooth solutions using a sequence 
of analytic solutions to a sequence of analytic "approximate equations." Although this 
method has successfully been applied by Clausen [20] to the family of polarized T 2 - 
symmetric spacetimes, and has also been used by Stahl [37] in studying S 3 and S 2 x S 1 
Gowdy spacetimes, it has proved difficult to apply in more general cases, such as for 
spacetimes with only one Killing vector field [26], or with no symmetries [4]. 

Our goal in the present work is to develop a general Fuchsian formulation for an- 
alyzing smooth (but not analytic) solutions to quasilinear hyperbolic PDEs which can 
be fairly directly applied to polarized and half-polarized T 2 -symmetric solutions to the 
Einstein equations and can be applied to polarized and half-polarized [/(l)-symmetric 
solutions as well. Two of the authors of this work, Beyer and LeFloch in [13], have 
carried out this program for semilinear hyperbolic systems and have applied their for- 
mulation to T 3 Gowdy solutions. Therein, Beyer and LeFloch set up a second-order 
Fuchsian formulation for smooth semilinear systems. In the present paper, in addition 
to generalizing to smooth quasilinear hyperbolic systems, we also work with Fuchsian 

In addition to the results we cite here involving quasilinear Fuchsian PDEs with non-analytic reg- 
ularity, there are a number of treatments of linear Fuchsian PDEs with non-analytic regularity; these 
include the works of Tahara [39, 38] and Kichenassamy [31, 28]. 
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systems in first-order form. 

One of the motivations for the semilinear work and its application to the T 3 Gowdy 
spacetimes was that the approximation scheme which plays a key role in the existence 
proof can also be used as the basis for a robust method for numerical simulations. 
This numerical approach has been developed and implemented in [12, 15] (see also [3]) 
as a tool for the numerical exploration of Gowdy solutions. Since our analysis in the 
present paper involves a similar approximation scheme, we expect to be able to carry 
out numerical investigations of singular initial value problems in more general classes of 
equations in future work. 

An outline of this paper is as follows. We begin the discussion of our results in 
Subsection 2.2, where we consider a general class of first-order quasilinear Fuchsian 
systems and, then, formulate the singular initial value problem for such systems. Next, 
in the same subsection, we introduce the class of first-order quasilinear Fuchsian systems 
in symmetric hyperbolic form and, for such systems, state an existence and uniqueness 
result (in Theorem 2.4, below) which holds in, both, infinite differentiability and finite 
differentiability classes. In Subsection 2.3, we carry out the details of the proof of this 
result. 

This theorem holds for a broad class of asymptotic data specified on the singularity. 
Next, in Subsection 2.4, we discuss special choices for this data, which we call "ODE 
leading-order term" , and then state and establish an alternative existence and uniqueness 
result (Theorem 2.21) for the singular initial value problem. This alternative theorem 
is useful in applications, as we show later in Subsections 2.4.4 and 3.3.3. Interestingly 
enough, ODE-leading-order asymptotic data also play a useful role as approximate 
solutions. 

In the second part of this paper, in Section 3, we apply our theoretical results and 
study polarized and half-polarized T 2 -symmetric solutions of Einstein's vacuum field 
equations. To this end, in Subsection 3.2, we define this family of spacetimes, and 
then, in the same subsection, we write out the Einstein equations in terms of areal 
coordinates. Next, in Subsection 3.3.1, we discuss the concept of AVTD behavior and 
discuss what one needs to do in order to check whether a set of polarized or half-polarized 
T 2 -symmetric solutions do exhibit AVTD behavior. In Subsections 3.3.2 and 3.3.3, we 
are in a position to rely on our results in Section 2 and we establish that our conditions 
therein indeed hold true for the spacetimes under consideration. On one hand, we rely on 
Theorem 2.4 and establish the existence of a parametrized family of polarized and half- 
polarized T 2 -symmetric solutions with, both, finite or infinite order of differentiability 
and with the expected AVTD behavior. On the other hand, we rely on Theorem 2.21 and 
show that, provided attention is restricted to smooth (C°°) solutions, then the family 
of T 2 -symmetric with AVTD behavior can be extended to include those with a wider 
("optimal") range of the "asymptotic velocities" (labeled by k, below, as defined and 
discussed in Section 3.3). 

In Section 4, we conclude and discuss the relevance of Fuchsian formulations for 
numerical simulations, as well as the application of the proposed formulation in order to 
tackle more general families of spacetimes. 



4 



2 The singular initial value problem 



2.1 Objective of this section 

For a given (say, first-order) PDE system V[ip] = 0, the (regular) Cauchy problem 
involves finding a solution tp = tp(t, x) to this system such that, at some chosen value to 
of the time, the solution satisfies the initial condition ip(t, xq) = <p(x), for some specified 
initial data function <p = <p(x). If the Cauchy problem is well-posed, then for any 
appropriate choice of (j), the solution ip exists for some open interval I containing to- 

Here, we are interested in the singular initial value problem rather than the regular 
one. That is, rather than seeking solutions to = which agree with specified 

initial data at a chosen time, we seek for solutions which become (in general) singular 
as t approaches some fixed value too, and which agree with some specified fall-off data 
as one approaches the singularity at too. As for the Cauchy problem, one can prove 
existence and uniqueness theorems for the singular initial value problem; these theorems 
guarantee that for any appropriate choice of the "asymptotic data" , there is a solution 
to V[tp] =0 which exists for t approaching too, and which asymptotically matches the 
prescribed asymptotic data 

2.2 Quasilinear first-order symmetric hyperbolic Fuchsian systems 

Before making this notion of singular initial value problem precise, and proving existence 
and uniqueness results, we carefully define the class of PDEs we shall consider here. 
While the theory we develop here can be generalized to a much wider class of background 
spacetime manifolds 1 (see, for example [2], in which we generalize our results to spacetime 
manifolds (0,5] x T n ), let us presume for now that we work on the cylinder spacetime 
(0,6] x T , for some small parameter 5, with coordinates t £ (0,5] and x £ T . The 
singularity is presumed to occur at t = (hence we set too = in the earlier notation); 
correspondingly, it is useful to work with the "singular time differential operator" 

D := td t . 

The general form of the first-order PDE systems under consideration reads 

S±(t, x, u)Du(t, x) + S2(t, x, u)td x u(t, x) + N(t, x, u)u(t, x) = f(t, x, u), (2.1) 

in which the unknown is a vector-valued spacetime function u : (0, 5] x T 1 — > R d for 
some integer d > 1 and some real 5 > 0. Here 5i = Si(t,x,u), S2 = S2(t,x,u), and 
N = N(t, x, u) are specified d x d-matrix~valued maps of the spacetime coordinates 
(t,x) and the unknown u (but is independent of its derivatives), while / = f(t,x,u) 
is a specified Revalued map of (t, x) and the unknown u (but again is independent of 
its derivatives). The specific requirements for the functions Si, S2, N and / are fixed 

1 We work with one spatial dimension, primarily since in addition to simplifying the discussion, this 
case is sufficient for handling our primary application: 3 + f dimensional spacetimes with a spatially 
acting T 2 isometry group. 
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precisely below; see in particular Definition 2.2. For notational convenience, we often 
leave out the arguments (t,x), instead we use the short-hand notation Si(u), 82(11), 
N(u), and f(u). Notationally, an object such as S\(u) may be interpreted as a map 
u 1 y Si(u) between two function spaces (as further discussed below). In this context we 
often write Si(u)(t, x), and we do the same for ^(u) and N(u). 

Observe that, in principle, one could absorb the term N(u)u into the source f(u); 
however in view of the conditions on these terms that we will introduce below, it is im- 
portant to keep these two terms separate. A system of this form (2.1) (noting especially 
the use of the singular operator D) will be referred to as a quasilinear first— order 
Fuchsian system. 

Before formulating the singular initial value problem for such systems, we wish to 
define functional norms and the corresponding function spaces which include built-in 
specifications of the asymptotic behavior of the functions in time. We state these defi- 
nitions first (here) for vector-valued functions, and then (below) for matrix-valued func- 
tions. The definitions are parametric: For the vector-valued function case, we specify as 
parameters i) a non- negative integer q, and ii) a fixed smooth 1 vector- valued function 
fi : T 1 — > M. d . Then using [i (which we label as an exponent vector) to construct the 
corresponding diagonal matrix 

n\ji](t,x) :=diag(r^*>,...,t-w(*>), (2.2) 

we define the norm 

IMk/x.q == sup \\K{fi](t,-)w{t,-)\\ Hq(T1) (2.3) 
te(o,<5] 

for vector-valued functions w = w(t,x). Here, || • H/j^t 1 ) denotes the standard g-order 
Sobolev norm on T 1 . Based on (2.3), we define the Banach space Xs^qiT 1 ) — also simply 
written as X§^ q — as the completion of the set of all functions w G C°°((0, 5] x T 1 ) for 
which this norm is finite, and we denote by Bg„q,r C X§n q the closed ball of radius 
r > (measured with the given norm) and center 0. To handle the class of functions 
which are infinitely differentiable, we define the space 

00 

In order to compare two function spaces Xs ifJ , t q and X^^,qi we write v > ji if, for each 
index i = 1, . . . , d and for all x G T 1 , the components of v and fi satisfy the inequality 
Vi(x) > im(x). Clearly, X SjV>q C X 8 ^ q if v > fi. 

We use analogously defined norms and functions spaces in order to control d x d 
matrix- valued functions such as Si, 52, and N. More specifically, in this case we choose 
as a parameter a fixed d x d matrix-valued valued function £ (labeled as an exponent 

x One could choose fi (and other parameter functions below) to have less regularity, but this is an 
inconsequential generality for the purpose of this paper. 
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matrix) which depends smoothly x £ T , and we define the corresponding norm as 

d 

\\S\\s, C , q := sup TTWt-teUSijiMWm, (2.4) 
te(o,s] fT± x 

for matrix- valued functions S = S(t,x). We denote the corresponding Banach space by 
Xg£ t q. Based on these function spaces of matrix-valued functions, we define (for r > 0) 
Bs,(,q,r C an d ^5,c,oo (as above), and we note that Xs^, q C Xg^ q if £ > C- 

As noted above, the singular initial value problem associated with a system such as 
Eq. (2.1) consists of choosing a set of "asymptotic data", and seeking for solutions which 
asymptotically approach that data. Using the function spaces just defined, we make this 
idea precise as follows. 

Definition 2.1. Given the parameters 6,fj, and q as above, and a chosen function uq : 
(0, 6] x T 1 — > R d , the singular initial value problem consists of seeking a solution 
u = uo + w to Eq. (2.1) whose remainder w belongs to Xg uq(T ). 

The function uq, which we refer to as the leading-order term, constitutes the 
asymptotic data, and is (a priori) of unspecified regularity. Regarding the solution 
function w, if it is to be considered a "remainder", then in comparison with uq it should 
be of higher order in t as one approaches t = 0. Observe that the exponent vector u, 
which parametrizes the function space Xg^^T 1 ), controls the order of the singularity 
of the remainder w at t = 0; roughly speaking, each component of w is of corresponding 
component order 0(V l ) if w G Xgn q . Hence, the components of \jl are sometimes referred 
to as the remainder exponents, with (i collectively labeled the (remainder) exponent 
vector. Generally, we assume here and below that exponent vectors are smooth. Also, 
for a given exponent vector a and a given scalar e, we use the notation a + e to indicate 
a new exponent vector obtained by adding e to each component of fi. 

We now discuss the conditions on Si(u), S^C"), N(u), and f(u) in (2.1) which, 
together with further conditions on the space of leading-order terms and the space of 
remainder functions, are sufficient for establishing the well-posedness of the singular 
initial value problem. The main set of conditions needed is included in the following 
definition. 

Definition 2.2. Fix positive constants 5 and s, a pair of non-negative integers qo and 
q (possibly +cc), and an exponent vector u : T 1 — > M. d , together with a leading-order 
term uq : (0,5] x T 1 — > M. d (with so far unspecified regularity). The system Eq. (2.1) 
is called a quasilinear symmetric hyperbolic Fuchsian system around uq if, for 
each x £ T 1 , there exist a matrix Si t o(uo)(x) that is positive definite and symmetric 
and independent of t, a matrix S2fi( u o)( x ) that is symmetric and independent of t, 
and a matrix Nq(uq)(x) that is independent oft, all defining matrix-valued functions in 
the Sobolev space H qo (T l ); and if there exists a smooth vector function (3 : T 1 — > M. d 
with strictly positive components, such that for every 5' G (0,(5], each of the "remainder 
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matrices " 

S\,i{uq + w) := Si(u + w) - S lfi (u ), 
5 2 ,i(u + to) := TZ\J3 - l]S 2 (u + to) - S 2i o(u ), 
JVi(«o + w) := N(u + to) - N {u ), 

considered as an operator of the form (for example) w t- > Si t i(uo+w), maps all functions 
w ^ Bs',u,q,s t° elements in Bgi £ q r} in which ( is some exponent matrix with strictly 
positive entries, and r > is some constant. It is furthermore required that 5i,i(tto + to) 
and 5*2,1 («o + w) are symmetric matrices for all to € Bg /xq,s- 

Before discussing further conditions which are needed in order to obtain existence 
and uniqueness for the singular initial value problem, we note the following remarks: 

(i) If a system Eq. (2.1) satisfies the conditions in the above definition and is thus a 
quasilinear symmetric hyperbolic Fuchsian system, then the matrices Si (it), S 2 (u), 
and N(u) in Eq. (2.1) (acting on u = uq + w) decompose as 

Si(u + w) = S lfi (u ) + Si^ituo + to), 

S 2 (u + w) = TZ[l - /3](S 2)0 (u ) + S 2j i(u + w)), 

N(u + w) = N (u ) + JVi(«o + w). 

Quasilinear symmetric hyperbolic Fuchsian systems are therefore "essentially" 
semilinear (described by the coefficients S^o, 5*2,0 and No), up to "quasilinear 
perturbations" (given by 5i,i, ^2,1 and N\), which decay as t — > with a rate 
controlled by £. The purely semilinear case has been treated earlier within a 
second-order framework by Beyer and LeFloch [12, 13, 14]. (See also [11, 15]). 

(ii) Presuming that a fixed leading-order term uo has been chosen, it is convenient 
to use the short-hand notation Si t \(w), S 2t i(w), Ni(w) in place of £1,1 (no + w), 
S 2 ,i(uo + w)i N\(uq + to), respectively, and similarly to use Si,o, 5*2,0, ^0 m place 
of 5i,o(tio), 5*2,0(^0)5 -^o(^o)- It is important, however, to keep in mind the depen- 
dence of these matrices on the choice of the leading-order term no . 

(Hi) If Eq. (2.1) satisfies the conditions in Definition 2.2, then it is symmetric hyper- 
bolic for all t S (0,5], provided 5 is sufficiently small in order to guarantee that 
5i is positive definite. Consequently, standard theorems ensure that the Cauchy 
problem with initial data specified at t = to G (0, 5] is well-posed in the usual sense 
(away from t = 0), so long as the order of differentiability (determined by q,qo) is 
sufficiently large. The solutions belong to the space C(J, ^(T 1 )) for q > 2 and for 
some interval I C (0, 5]; however, nothing is known a priori regarding the behavior 
of these solutions as t approaches the singularity at t = 0. 

(iv) We note that the matrix-valued operators 51,1(10), 52,i(to), and N\(w) are non- 
singular in a neighborhood of t = since they take values in Xg^ iq with £ > 0. 
Hence quasilinear symmetric hyperbolic Fuchsian systems are singular precisely at 
t = 0, where the PDE coefficients are singular. 
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(v) Presuming that q is sufficiently large, the conditions in Definition 2.2 require that 
the coefficient matrices on the left-hand side of Eq. (2.1), e.g., S\(t,x,u), are 
defined on the domain (0, 5] x T 1 x U, where U is an open subset of M. d about the 
origin. This domain U must be compatible with the choice of parameters 6, [i and 
s. For example, if fj, > 0, q > 1 and w G Bg jlJi>qjS , then < Cs5 Mmin , where 

C is the Sobolev constant and // m i n is the minimal value over all components of // 
over all spatial points. Hence, if necessary, s and or 5 must be chosen sufficiently 
small in order to fit into U. 

We discuss a collection of useful technical tools in the appendix (primarily in Sec- 
tion B) which allow us to check if the conditions of Definition 2.2 are satisfied for a given 
problem. 

The remaining conditions we consider concern the coupling between the components 
of an valued function u = u(t,x), presumed to satisfy Eq. (2.1), and the effects 
of these couplings on the asymptotic behavior of the components as t approaches the 
singularity. 

Definition 2.3. Given the singular initial value problem (Definition 2.1) for a speci- 
fied quasilinear symmetric hyperbolic system (Definition 2.2) with specified leading-order 
term uq and specified function space Xgn q , the system Eq. (2.1) is called block diagonal 
with respect to fi, provided the following commutation conditions hold 

IZ^S^u) = Si(«)ft[/x], K[fj]S 2 (u) = S 2 (u)K[p], lZ[y]N(u) = N(u)K[fi}, 

for all u = uq+w with w G Xs ^^ q , and provided the same condition holds for all relevant 
spatial derivatives of Si(u), S2(u), and N(u). (Recall that TZ[fi] is defined in Eq. (2.2).) 

This block diagonality condition is used in the derivation of energy estimates (in 
Section 2.3, below) and thus plays a major role in the proof of existence and uniqueness. 
Roughly speaking, this condition guarantees that the "principal part operator" 1 

L(u)[v] := S!(u)Dv + S 2 (u)td x v + N(u)v, (2.5) 

takes block diagonal form for v, and that each block is associated with only one com- 
ponent of the exponent vector /i. Recall that the components of \x determine the order 
of the singularity in t at t = for the components of the remainder of the singular 
initial value problem. Hence, this condition requires that the principal part may only 
be coupled within those components of the solution whose remainders behave the same 
at t = 0. Note that the condition does allow all of the matrices S\, S 2 and N to depend 
on all components of u. It also allows for arbitrary coupling in the source-term, which 
we write as 

F(u )[w] := f(u +w), 
(or, in short form, as whenever it does not lead to confusion). 

1 This operator appears in Eq. (2.1) in the form L(u)[u]; however, it is useful in the discussion below 
to define this operator with distinct functions u, v. 
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We are now ready to state our main existence and uniqueness results for the sin- 
gular initial value problem associated with first-order quasilinear symmetric hyperbolic 
Fuchsian systems. Our hypotheses below include conditions on the matrix 

M := Si i0 DiagOui, . . . , fj, d ) + N , (2.6) 

which we refer to as the energy dissipation matrix and depends on the space coor- 
dinate x, only. 

Theorem 2.4 (Existence theory for symmetric hyperbolic Fuchsian systems). Suppose 
that Eq. (2.1) is a quasilinear symmetric hyperbolic Fuchsian system around a leading- 
order term uq (with the choice of the parameters 1 5, s, \x, q, and qo as specified in 
Definition 2.2) and is block diagonal with respect to fj,. Suppose that q > 3 and qo = q+2. 
Then there exists a unique solution u to Eq. (2.1) whose remainder w := u — uq belongs to 
with Dw G ^ j for some 5 £ {0,5], provided the following structural conditions 
are satisfied: 

(i) The energy dissipation matrix Mq defined in Eq. (2.6) is positive definite for every 
spatial point iGT 1 . 

(ii) The map 

■J{u G ) : w i— > F{u ) [w] - L{u + w) [no] (2.7) 

is well-defined, and for every 5' G {0,5], it maps w G Bgi^ qjS to X$i uq for some 
exponent vector v > \i. 

(Hi) For all 5' G (0, 5], for a constant C > and for all w,w G Bgr u q s , one has 

|| JH - 3[w] \\ s , iVA <C\\w- w\\ff, M (2.8) 

and 

\\3[w] - ^[w]\\ S ',u,q-l + \\Sl,l{w) - Sl,l(w)\\s> £, q -l 

+ \\S2,l{w) ~ S 2 ,l{w)\\s',C,q-l + \\Nl{w) - Nl(w)\\s' £, q -l < C\\W - w\\s> 

(2.9) 

for all w,w G B S ',^q, 8 - 

If all of these conditions are satisfied for all q > 3, then there exists a unique solution u 
of Eq. (2.1) such that u — uq and D{u — Uq), both, belong to 2 ^j MOO - 

1 Although Definition 2.2 involves the choice of j3, ( in addition to 5,/j,,q,qo, this latter set is crucial 
to the nature of the results of this theorem, while the former set is in a certain sense incidental. Hence 
even though £ does appear in condition (iii), we do not list /3 and £ in the hypothesis here. 

2 This result is short of "well-posedness" in the conventional sense because we do not prove continuous 
dependence of the solution on the asymptotic data. Such a result is certainly desirable, and we plan to 
investigate it in future work. 
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Section 2.3, below, is devoted to the proof of this theorem. Observe that, in the 
hypothesis of this theorem, the regularity required for S^o, S^o, and Nq (specified by 
qo) slightly differs from that required for S\^, S^i, and N\ (specified by q). The same 
observation can be made regarding the asymptotic data uq (implicitly specified by Con- 
dition (ii)) and the solution u (specified by q). These gaps arise in the course of our proof, 
in particular in obtaining the energy estimates for the Cauchy problem Lemma 2.8. It is 
not clear whether this discrepancy in regularity could be eliminated by another method 
of proof, and in any case it disappears in the "smooth" case, if q and go are both infinite. 

In formulating this theorem, we require that the source term operator w t— > 3~[u>], 
and with it the source term function f(t,x,u) in Eq. (2.1), be defined on the domain 
(0, 6] x T 1 x U, where U is an open neighborhood of the origin in M, d . In the same way as 
for the coefficient matrices Si, S2 and N, we find that the parameters 5, \x and s must 
be compatible with U. We note that Condition (ii) also restricts the regularity of the 
leading-order term uq. 

We also note that the time of existence of the solutions, specified by 5, could a priori 
be very small. Indeed, a smaller choice of the parameter s (which may be necessary in 
order to fit into the domains of the coefficient functions of Eq. (2.1)) generally leads to 
a shorter guaranteed time interval of existence. 

In its applications, Theorem 2.4 often allows one to find an open set of values for the 
exponent vector fj, for which the singular initial value problem admits unique solutions. 
A lower bound for this set 1 can originate in Condition (i), while an upper bound is 
usually determined by Condition (ii). Both bounds on the set of allowed values for fi 
provide useful information on the problem. The upper bound for /i specifies the smallest 
regularity space and, hence, the most precise description of the behavior of w (in the limit 
t \ 0)i while the lower bound for \i determines the largest space in which the solution u 
is guaranteed to be unique. Observe that this uniqueness property must be interpreted 
with care: under the conditions of our theorem, there is exactly one solution w in the 
space Xc , although we do not exclude the possibility that another solution may exist 
in a larger space, for example, in ~ with ju < \i. Note that if a given system does 
not satisfy our hypothesis above, there is sometimes a systematic method which allows 
one to "improve" a leading-order term uo; cf. the discussion of (order-n)-leading-order 
terms in Section 2.4 and, in particular, Theorem 2.21. 

We also remark that results analogous to those stated in Theorem 2.4 for D(u — uq) 
can be derived for higher-order time derivatives of the solution. 

2.3 Proof of the existence and uniqueness theorem 
2.3.1 Outline of the argument 

Before carrying out the details of the proof of the existence and uniqueness Theorem 2.4, 
we outline the basic strategy and the basic steps of the proof. We start by working with a 

A real A is defined to be a lower bound for the allowed values of the vector if each component v a 
of v satisfies the condition v a (t,x) > A for all x in the domain of fi. A similar definition holds for an 
upper bound for [i. 
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linear version of the PDE system. We consider the Cauchy problem for this linear system, 
verifying that the conditions we have assumed as part of the hypothesis of Theorem 2.4 
guarantee local existence and uniqueness of solutions for this Cauchy problem, with 
appropriate levels of regularity. We then use these results pertaining to the Cauchy 
problem for the linear system and establish that unique solutions to the singular initial 
value problem for the linear system exist in a neighborhood of the singularity. This is 
done using the solutions of sequences of Cauchy problems with the initial time for the 
j'th element of this sequence set at tj, and with tj approaching zero, the time of the 
singularity. To show that the limit of such a sequence of solutions exists, and satisfies 
the singular initial value problem, we work with the linear PDE system in a weak form, 
and we also employ a family of energy functionals. To proceed from solutions of the 
singular initial value problem for the linear system to solutions for the full quasilinear 
system of Theorem 2.4, we use a standard fixed point iteration argument for a sequence 
of linearized equations and their singular initial value problems. Observe that arguments 
similar to those used here have been applied in [12] in order to establish existence and 
uniqueness results for the singular initial value problem for semilinear (second order) 
Puchsian PDEs. 

2.3.2 The singular initial value problem for linear PDEs 

The linear systems we consider here are essentially those of Theorem 2.4 (see Eq. (2.1)) 
with Si, 52, and N set to be independent of u, and with / set to be linear in u. More 
specifically, we introduce the following definition. 

Definition 2.5. Suppose that 5 and r are positive reals, q and qo are non-negative 
integers, fi : T 1 — > M. d is an exponent vector, and £ : T 1 — > M. dxd is an exponent matrix 
such that ( > 0. The system Eq. (2.1) is called a linear symmetric hyperbolic 
Fuchsian system if the following conditions are satisfied: 

(i) The operators Si, £2, and N are independent of u, and they can be decomposed as 

S 1 (t,x) = S 1>0 (x) + S 1>1 {t,x), (2.10) 
S 2 (t,x)=Tl[l-P(x)]{s2,o(x) + S2 t i(t,x)), (2.11) 
N(t,x) = N (x)+N 1 (t,x), (2.12) 

where Si$ is symmetric and positive definite at every spatial point, where S± t i, S2,o, 
and 62,1 are symmetric, and in addition, the maps Si 5 o,S2,o and Nq belong to 
H qo (T ), while Si,i,S2,i and Ni are d x d-matrix-valued functions in Bs^, q , r - 
Here, ft : T 1 — > R rf is a smooth vector function with strictly positive components. 

(ii) The constant 5 is sufficiently small so that Si is uniformly positive definite 1 . 

1 By uniformly positive definite we mean that Si(t,x) = 5*1,0(2;) + Si,i(t,x) is uniformly positive 
definite (I) at all t E (0,5] in the L 2 -sense (with respect to x), and (11) for all Si,i E Bsx,q,r- Since 
we assume in addition above that Si t o is positive definite and the perturbation Si,i is bounded, this 
implies that Si is positive definite pointwise at all (t, x) E (0,(5] x T 1 for all Sij £ Bsx,g,r if 9o and q 
are sufficiently large (thanks to the Sobolev inequalities). 
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(Hi) The source term is linear in the sense that 

F(u )[w] = f(t,x,u + w) = f (t,x) + Fi(t,x)w, 

with /o G ^S,v,q and the matrix F\ satisfying 1Z[ll]FiR[[i\~ 1 G 13$^^^. Here v is 
an exponent vector with v > 

In this definition, we note the condition v > [i. It is used primarily in the proof of 
Proposition 2.10, to enforce the needed rapid decay of the source term fo(t, x) as t — > 0. 

It is clear from this definition that a linear symmetric hyperbolic Fuchsian system is 
a special case of a quasilinear symmetric hyperbolic system, with the leading-order term 
uo = (this is no loss of generality for linear systems). Both in the linear and in the 
quasilinear case, we consider the functions S^i, $2,1, N±, and F\ to be perturbations of 
S\fi, 5*2.0, -^0) an d fo- An important step in our analysis is to seek uniform estimates for 
these perturbations. It turns out that such estimates can only be obtained if the pertur- 
bations are bounded. This is the reason for introducing the balls with radius r above, 
Bs£,qri which can be considered as those spaces in which we seek the perturbations. 

In carrying out the proof, it is important that we keep careful track of which quan- 
tities the constants arising in various estimates are allowed to depend upon. To make 
this precise, it is useful to have the following definition. 

Definition 2.6. Suppose that Eq. (2.1) is a linear symmetric hyperbolic Fuchsian system 
for a chosen set of the parameters 5, u, £, q, qo and r. Suppose that a particular estimate 
(e.g., the energy estimate Eq. (2.16)/, involving a collection C of constants, holds for 
solutions of Eq. (2.1) under a certain collection of hypotheses H. The constants C are 
defined to be uniform with respect to the system and the estimate so long as the following 
conditions hold: 

(i) For any choice of Si ; i, S^i, Ni and F\ contained in the perturbation space Bg^g r 
(see Definition 2.5) which is compatible with the hypothesis H, the estimate holds 
for the same set of constants C. 

(ii) If the estimate holds for a choice of the constants C for one particular choice of 
5, then for every smaller (positive) choice of 5, the estimate remains true for the 
same choice of C. 

Recalling our definition above (see Eq. (2.5)) of the principal part operator L, we 
define the linear principal part operator by 

L[w] := (Si,o + S!,i)Dw + n[l - P](S 2 ,o + S 2 ,i)td x w + (N + Njw. (2.13) 

In terms of this operator, the linearized PDE system Eq. (2.1) may be written in the 
form L[w] = fo + F\w. 

In summary, the parameters 5, \x and q determine the space Xg^q for the remainder 
of the solution of the singular initial value problem with leading-order term uo, while S, 
C, q and r fix the space Bg^ q>r of the perturbations of the coefficients. The parameter go 
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determines the order of differentiability of the "leading-order" coefficient matrices 5i,o, 
5 2 ,o and No- 

Suppose that Eq. (2.1) is a linear symmetric hyperbolic Fuchsian system (for param- 
eters 5, r, q, qo, £, fi; cf. Definition 2.5). We first consider the Cauchy problem; that 
is, we prescribe initial data vu ] specified at some to G (0, d~) and we seek solutions on 
[to, 6] x T 1 which agree with vu Q -\ at t = to- It is useful at this stage for us to make 
the temporary assumption that 5i,i, 5*2,1, -^l an d F\ are C°°((0, 5] x T 1 ) functions con- 
tained in their respective function spaces, as discussed in Definition 2.5. (If Si i, 5*2,1, 
Ni and F\ satisfy this smoothness assumption, then Eq. (2.1) is said to have smooth 
coefficients 1 .) Given such a linear symmetric hyperbolic system Eq. (2.1) with smooth 
coefficients and if in addition qo > 2 and also fo G ^s,v,q is smooth, then it is a standard 
result (see, e.g., Chapter 16 in [40]) that the Cauchy problem is well-posed in the sense 
that for initial data vu ] G H q °(T l ), there is a unique solution v : [to, 5] xT 1 — > R rf to 
this Cauchy problem with v(to) = Vu Q i and with v(t, •) G H qo (T l ) for all t G [to,£]- It 
is crucial for the following discussion that indeed this solution exists for the full inter- 
val [to,$], regardless of the choice of to G (0,(5). This is true as a consequence of the 
positivity of Si on (0,5]; cf., Condition (ii) in Definition 2.5. 

In fact, this statement about the Cauchy problem remains true if the matrices 5i.i, 
52,i, N±, fo and F\ are not required to be smooth, but are only required to have qo spatial 
derivatives. Such a relaxation is, however, not useful for our arguments; we use a more 
general continuation argument below by which we recover the non-smooth case. We also 
note that although this assumption that Eq. (2.1) has smooth coefficients implies that 
5i,i, 52,i, N\, fo, Fi are differentiable to all orders, it does not not guarantee that all 
derivatives have controlled asymptotic behavior for t \ 0. This control holds only for a 
set of derivatives given by q, as labeled by the relevant function space. 

To establish control over the solutions to the Cauchy problem for the linear version of 
the PDE (2.1) and the regularity of these solutions, we now introduce a two-parameter 
family of explicitly time-dependent energies: Presuming that the remainder exponent 
vector \x is fixed, for any pair of positive constants k and 7 we define the energy Ea K ,y 
for a function w : [to, 5] x T 1 — >■ R d (with w(t, •) G L 2 (T X ) for each t G [to, 5]) as follows: 

E^AMit) ■= \e~ Ktl (5i(t,-)^[A*](t,-Mt,-),W(t, •)«;(*, 0>£a(Ti)> ( 2 - 14 ) 

where 5i is the matrix appearing in Eq. (2.1). We emphasize again that, unlike standard 
definitions of energy, the energy functionals Ea yK ~\w\(t) defined here depend on time ex- 
plicitly, and not just through the time dependence of w(t, x). Note that it readily follows 
from this definition, and from the conditions assumed to hold for 5i in Definition 2.5, 
that there exist uniform (in the sense of Definition 2.6) constants C\ and C2 such that 
for any L 2 {T l ) function w{t,x), one has (for all t) 

CxWTl^it, >(t, OHl^) < ^E^[w](t) < C 2 \\n\p\{t, >(t, Oliver*)- (2-15) 

1 Since the assumption of smooth coefficients is not presumed to hold everywhere below, in all cases 
that is does hold, we state that explicitly. 
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These energies, as is usually the case, have been defined in such a way (including 
the presence of the factor e~ KtJ ) that for solutions of the Cauchy problem for Eq. (2.1), 
which we label v(t,x), the growth of the energies is controlled. Explicitly, we obtain the 
following estimate. 

Lemma 2.7 (Basic energy estimates for the Cauchy initial value problem). Suppose 
that for some choice of the parameters 5, fi, (,q,qo, and r, with q = and qo = 2, and 
for uq = 0, Eq. (2.1) is a linear symmetric hyperbolic Fuchsian system with smooth 
coefficients and with /o both smooth and contained in X$^ q for some v > fi. Suppose 
also that the system is block diagonal with respect to n, that the energy dissipation matrix 
Eq. (2.6) is positive definite for all x £ T 1 and, in addition, that DSi t i and d x S2i are 
contained in Bggog for some constant s > and some exponent matrix £ with strictly 
positive entries. Then there exist positive constants k, 7, and C such that for any initial 
data i^ ] E ff 2 (T 1 ) specified at some to G (0,5], the solution of the Cauchy problem v 
for this system and this initial data satisfies the energy estimate 

y/E^[v]{t) < \J E^ K>1 [v\(t)\ t = t0 +C s~ 1 \\11[li]{s, -)/o(s, Ollz^CT 1 )^ 5 (2-16) 

J to 

for all t £ [fo)£]- The constants C, k, and 7 may be chosen to be uniform 1 and do not 
depend on f$. In particular, if one replaces vu^ specified at to by any ur tl i specified at 
any time ti G (0, to]) then the energy estimate holds for the same constants C, k, 7. 

Before proving this lemma, we make a few remarks: I) Lemma 2.7 does not imply 
that the energy estimate Eq. (2.16) holds for t < to; in particular, it need not hold 
for t \ 0. II) The well-posedness of the Cauchy problem which is used implicitly in 
the proof of Lemma 2.7 requires sufficiently high regularity on the coefficients (see for 
example [40]); this gives rise to the condition qo = 2 stated in the hypotheses. Ill) We 
remind the reader that the condition that the coefficients be smooth does not imply 
either the q = condition or the conditions that DS\ t \,d x S2,\ £ -65,^,0,5- While the 
smoothness condition implies the existence of all derivatives, the latter are statements 
about the behavior of the lowest derivatives in the limit t \ 0. It may appear that since 
Lemma 2.7 focuses on the Cauchy problem at times to > only, control of behavior 
near t = is not necessary. However, such control is in fact needed to obtain an energy 
estimate with constants which are uniform and independent of to- IV) In view of the 
norm equivalence (2.15) stated above, the estimate Eq. (2.16) can be rewritten as 

\\n[fi](t, >(t, Oilmen < c(\\K[n\(t , •Klb(Ti) 

ft l , (2-17) 

+ / s~ 1 \\n[^}(s r )fo(s,-)\\ L 2 {T1) d s ). 

We observe that for this version of the energy estimates, all of the constants C, n and 
7 are absorbed into the constant C; every change of the former constants is therefore 

1 While the constants C, k and 7 here can be chosen to be uniform in the sense of Definition 2.6, 
there generally does not exist a choice which holds for all 5, Si,o, S2,o, No, /3, r, s, £, fi and v. 
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reflected in a corresponding change of the latter one. V) For some of the following 
discussion it is important to note that the particular values of the parameters of the 
perturbations space £ and r (and also £ and s) do not play an essential role in this 
lemma: if we change from one perturbation space Bg^^^ to another one B s ^ ~ (and £ 
and s), the same result is obtained with possibly different, but still uniform, constants 
C, 7, k. This is true for all of the following results. 

Proof. The basic idea of the proof is to compute DE[v](t), then bound the terms on the 
right hand side and finally integrate the equation in time. For simplicity we write E[v] 
in place of E^ iKi j[v]. Computing 1 DE[v], and using the symmetry of the matrix Si, we 
obtain 

DE[v](t) = - K^tf E[v]{t) + \e~ Ktl [ ((DSi)n[fj]v,n[fj]v) dx 

2 J T i 

+ e- KtJ [ (S^DK^v^K^v) dx + e~ RV [ (SiR\p\ Dv, K\p\v) dx. 
Jt 1 Jt 1 

We first analyze the fourth term on the right hand side of this expression, which we label 
/. Using the fact that v is a solution of equation Eq. (2.1), using the block diagonal 
condition (Definition 2.3), and integrating by parts, we calculate 

I = e~ Ktl J ( {n\ji]f ,n\ji]v) + {n[ti}Fiv,n[»]v) + l/2t {(d x S 2 )K\j\v,K[v\v) 

+ t (S 2 (d x Tl[fj])v,n[fj]v) - (Nn\jj]v,Tl[fj]v)^dx. 

Using the Holder inequality, we may then estimate the first term in this expression as 
follows: 

e~ Ktl [ (^[/i]/o,^[M]«)dx<e- Kt7 ||^M/ |U 2 ||^M^|U 2 . 
Jt 1 

We now argue that for appropriate choices of k and 7, all the other terms besides 
this one can be neglected in a certain sense. First, we use the properties of the linear 
symmetric hyperbolic Fuchsian system to expand the coefficient matrices Si, S 2 , N into 
terms which are O(l) at t — > 0, and terms which decay as a power of t. We thereby 
obtain 

DE[v] < - e~ KtJ [ ((N - SifiiDKin^Klfj]- 1 ) K[fj]v, K\jj]v) dx 
Jt 1 

+ e~ Ktl J ti ((-I k ^Si + S^DTZ^n^]- 1 + ±DS!,i - N x + K^F-frfr]- 1 

+ ^td x S2 + ts 2 {d x n[n})n[fi}^TZ[fi}v,n[fi}v^dx 

+ e- Ktl \\n[fi]f \\ L 2\\TZ[fi]v\\ L 2, 

In calculating this time derivative, we use the fact that the solution v is C 1 in both time and space. 



16 



where we use the expansion for 5 2 to write 

td x S 2 = d x n[-P(x)] (5-2,0 + S 2 ,i) + K[-p{x)] (d x S 2 ,o + d x S %1 ) . 

The first integral on the right hand side of this inequality for DE[v] is negative 
definite if the energy dissipation matrix Mq = Nq — Sifi(DTZ[/j])Tl[fi\~ 1 (see Eq. (2.6)) 
is positive definite, and so can be neglected. All of the terms in the second integral on 
the right hand side of this inequality decay as some positive power of t. We also know 
that as a consequence of Definition 2.5, the matrix S\ is positive definite uniformly. It 
is at this point that we use the factor of e~ Ktl which appears in the definition of the 
energy functionals. The scheme is to choose k and 7 in such a way that the second 
integral in the estimate above is negative definite. This can be achieved if we choose 
7 small enough and k large enough so that the negative definite Si-term dominates all 
of the other terms in the second integral on (0,5]. To see that the constants k and 7 
may be chosen so that they are independent of the functions 5i,i, £2,1, N\ and F\ and 
are therefore uniform in the sense of Definition 2.6, we recall that by assumption (see 
Definition 2.5), 51,1,52,1 an d N\, and TZ[^\FilZ[^\~ l are contained in the ball Bs^ iq>r . 
Hence these functions all must have finite norms bounded by r. Since the role played by 
5i,i, 52,i, N\ and F\ in determining the constants k and 7 depends strictly on the norms 
of these functions, we may choose a fixed set of the constants such that the inequality 
holds for any 5i,i, 52,i, N\ and F\ contained in these balls. In total, we obtain 

de[v] (t) < e - Ktl 1 \niiAM\u2 1 \n[fi]v\\ L 2 , 

which implies that 

d t E[v](t) <r l e-^\\Tl[»}M\ L 2\\n[»]v\\ L2 . 

Then using the norm equivalence Eq. (2.15), we may rewrite this as 

d t E[v](t) ^Ct-^-^II^M/oll^v^M- (2.18) 

To integrate this inequality, it would be useful to divide both sides by y/E[v](t). 
However, since the L 2 norm of v may vanish in special cases, we use the following 
strategy. We set E e := E + e for some constant e > (see, for instance, [36, Page 
59]), and we check that (2.18) holds if we replace E by E e . Then dividing, and using 
-krd t E e = 2d t \fE~l, we obtain 

V -C'e 

d t ^/E^m <cr l e- Ktl \\nw Q \\ L 2, 

after rescaling the constant C . We now integrate both sides over ds, thereby obtaining 
VEM^)<VEe[v](to) + C /VV" s7 ||ft[//]/o|| L2 ( S )d S 

J t 

< y/E e [v](to) + c( sup e" Ks7 ) / s- 1 \\n[^]fo\\ L 2(s)ds 

v se(t ,t) 7 J to 

< ^/E € [v](t ) + C f s-'WKWoWMds, 

J to 
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where we note that the constant C changes from the second to the third line of this 
calculation. Taking the limit e — > finishes the proof that the inequality (2.16) holds. 
It also follows directly that the constant C is uniform. □ 

In order to derive the solution of the singular initial value problem from a sequence 
of solutions of the Cauchy problem, we need estimates involving higher order spatial 
derivatives. We obtain these as follows. 

Lemma 2.8 (Higher-order energy estimates for the Cauchy initial value problem). Sup- 
pose that a linear symmetric hyperbolic Fuchsian system has been chosen which satisfies 
all of the conditions of the energy estimate Lemma 2.7, except that 1 (rather than q = 
and qQ = 2) q is an arbitrary integer greater than one, and qo = q + 2. Then there exists 
a pair of positive constants C and p such that for every sufficiently small e > 0, the 
solution v of the Cauchy initial value problem with initial data i^i € H qo (T l ) specified 
at to satisfies (for all t £ [to, 5]) 

\\n\fi - e](t, -)v(t, -)|k«(Ti) < C \\H[l4(to, >toll#9(Ti) 

rt (2 19) 

OII-Hs-ifr 1 )) ds. 

J t 

The constants C (which in general differs from C in Eq. (2.17)) and p are uniform? 
in the sense of Definition 2.6 and do not depend on /q. If we replace v\ to \ specified at 
to by any vu^ specified at any t\ G (0, to], then the same estimate holds, for the same 
constants C and p. 

Observe that it is necessary (as stated in the hypothesis of this lemma) that the 
solution (and hence the data and coefficients) be contained in H q+2 if we wish to obtain 
an energy estimate for q spatial derivatives. The reason for this requirement is made clear 
in the course of the proof. The main difference between the hypotheses of Lemma 2.7 
and 2.8 is that we require stronger control of the behavior of spatial derivatives of S\, 
$2, N, /o, and F\ in the limit t \ in Lemma 2.8 (i.e., q > 1 as opposed to q = in 
Lemma 2.7), while we presume smoothness for Si, S2, N, fo, and F\ in both lemmas. 

Proof. This lemma is proven by taking q spatial derivatives of Eq. (2.1), reorganizing the 
resulting equations into a linear symmetric hyperbolic Fuchsian system for the g'th order 
derivative of v , applying Lemma 2.7 to that system, and then carrying out a number of 
estimates needed to derive Eq. (2.19) from the inequality resulting from this application. 
We discuss some of the details for the q = 1 case here; the q > 1 cases are similar. 

We need the conditions d x S2,i £ Sa,j,o,s and DSi,i £ Bs,^,o, s to both be implicitly included in the 
hypothesis for this lemma. The first of these follows automatically from the following assumptions on 
q and qo- The second does not, but it is included in the hypothesis for Lemma 2.7, and therefore is 
included in the hypothesis for this lemma. 

2 The constant C generally only depends on 8, 5*1,0, 6*2,0, M), f3, r, fi, e and q, and the constant p 
depends only on e and q. 
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Presuming that v(t, x) is the solution to the Cauchy problem for the linear system 
(before differentiating) with initial data vi to ] (contained in H^lT 1 )) specified at to, we 
carry out the differentiation and obtain the following PDE system for d x v: 

SxDd x v + S 2 td x (d x v) + Nd x v = f + F x d x v, (2.20) 

where 

To ■= d x fo + {d x F x - d x N)v + d^S^iNv - /„ - F lV ), (2.21) 

and 

F, := Fx - ta x s 2 + ta x s 1 Si 1 s 2 . 

Here, we interpret v as a given function so that /o can be considered as a source term 
function, and d x v as the unknown. This PDE system is clearly of the desired form 
Eq. (2.1) (with d x v G H qo (T l ) for each value of t). However, it is not a linear 
symmetric hyperbolic Fuchsian system with respect to the same exponent vector fi: 
the term d x Nv in /n violates Condition (iii) of Definition 2.5 since it is in Xg^q-i rather 
than in Xg V q—\ for some v > \i. However, Eq. (2.20) is a linear symmetric hyperbolic 
Fuchsian system if we choose ]2 := fi— e/2 as the remainder exponent vector for any scalar 
constant e > 0. One verifies that Eq. (2.20) has block diagonal form with respect to /I 
and also that the energy dissipation matrix is positive definite if e is sufficiently small. 
Consequently, this system Eq. (2.20) satisfies the hypothesis of Lemma 2.7. It follows 
that there exist uniform (in the sense above) constants C, k and 7 (which generally differ 
from the ones for the original equation) such that d x v satisfies the energy estimate (for 
all t G [t , 5}) 

^E n ^[d x v](t) < J E % ^[d x v]\ t0 +C f s- l \\Kfr]{s, -)Us, -)|U a( Ti)<fe. (2-22) 

J to 

To derive the q = 1 version of the estimate (2.19) from the energy estimate (2.22), 
we first note two useful inequalities. Letting / : (0, 6] x T 1 — > M. d denote any function 
for which the following norms are finite, we find 1 that, for all t G (0,5], 

1 

\\K\fL - e]f(t)\\ HHTl) <CY, \\n\ji - e/2]4f(t)\\ L2{T iy (2.23) 

«=o 

Here, the constant C > may depend on \x and e, but, in particular, is independent of 
t. Observe that the use of fi — e on the left hand side of Eq. (2.23) and of \i — % on the 
right hand side, is needed to dominate the terms of the form log t which are picked up 
on the left hand side when R\p] is differentiated in space if \i is not constant (as a result 
of the ff 1 (T 1 ) norm). We also readily check that 

\R\ii - e/2]d x f(t)\\ LHTl) < C\\n[»]f(t)\\ H i iT iy (2.24) 

1 There are versions of this inequality, as well as the one below, which hold for higher order spatial 
derivatives of / as well. 
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We now work on inequality (2.22): Observe first that inequality (2.16) from Lemma 
2.7 holds for /2 = fj, — | so long as e is sufficiently small; hence we may add the left 
hand side of Eq. (2.16) (with fi — I) to that of Eq. (2.22), and the right hand side of 
Eq. (2.16) (again with /i— § ) to that of Eq. (2.22). If we now use (i) the norm equivalence 
Eq. (2.15) on both sides to replace energy terms by terms involving norms of 1Z[-], (ii) 
the definition of the Sobolev norm || • to combine terms on each side, and (iii) 

the inequalities (2.23) and (2.24), then we obtain the following inequality: 

||ft[//-eH H i(Ti)(t) <C^\n[fj]v\\ m{T i ) (t )+ £ S " 1 (||7e[ / u-e/2]( S ,-)/o(s,-)||L2(Ti) 

+ \\K[ti-e/2](s,-)f (s,-)\\ LHT t ) y s j. 

It remains to substitute in the definition of /o from Eq. (2.21). Noting the properties 
of the functions on the right hand side of Eq. (2.21), we determine that there exists a 
uniform constant C (in the sense above) such that, for all t £ (0,5], 

\\Kfr - e/2]/o(t)|| ia(T i) < C - e/2]d x f (t)\\L*(P) + \\K\m ~ e/2]v(t)\\ L2{Tl) ) 

L 2 (T 1 ) ) ■ 

Here in the second step, the constant C has been inconsequentially changed. Combining 
these last two inequalities, we obtain the desired result Eq. (2.19) with q = 1 by setting 
p = e/2. 

This concludes the proof that the inequality (2.19) holds for the case q = 1. The 
proof for q > 1 proceeds very similarly. The argument that the constants C and 5 may 
be chosen so that the inequality holds for all Si s i, 52,1, Ni, and Fi contained in Bs t c„q,r 
is essentially the same as that used in proving Lemma 2.7. □ 

We remark that while the introduction of e into the estimate Eq. (2.19) is certainly 
needed, one can choose this e to be arbitrarily small. One might worry that as one 
proceeds from q = 1 to higher values, the necessary value of e grows and causes trouble. 
However, since the incremental value needed for each step is arbitrarily small, one sees 
that the total value of e needed for arbitrary differentiability values can be kept small 
(below any chosen positive value). 

With these results for the Cauchy problem for linear symmetric hyperbolic Fuchsian 
systems established, we now set out to use solutions of the Cauchy problem to establish 
the existence of solutions to the singular initial value problem. We do this via an ap- 
proximation scheme which works as follows: We first choose a monotonically decreasing 
sequence of times t n £ (0, 5] which converges to zero. Then for each n, we construct a 
function v n : (0,5] x T 1 — > M. n which vanishes on (0, t n ], and which is equal on (t n ,5] 
to the solution of the Cauchy problem with zero initial data at t n . One readily checks 
that for every choice of \i, one has v n G C°((0,5] x T 1 ) n X s ^o- The central result of 
this section is that if certain hypotheses hold, then the sequence (v n ) - whose elements 
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we label approximate solutions - converges to a solution of the singular initial value 
problem for the linear system with vanishing leading term. 

The first step in showing this convergence is to set up the formalism to work with 
weak solutions to the linear system. To do this, we define a test function for this 
system to be any smooth function <fi : (0, S] x T 1 — > R d for which there is a J 1 £ (0, 5], 
such that 4>(t, x) = for all t > T. We then define the operators C and T acting on 
functions w E -^5,^,0 via 1 

(C[w],(j>) := - J (^(llinlSxw, D(p) L2{T1) + (TZ[fj]S 2 w,td x (p) L 2 [Tl) 

+ (n[fj] (Si - N W + n^Dni^S! + ds! + n[^}- l td x iz[^s 2 + td x s 2 ) , <t>) L2{T1) )dt 

and 

I' 8 

' l) L 2 {T 1 



rS 

(F[w], </>):= I (%](/ + Fit«)» i2(Tl) (it, 





where (p is an arbitrary test function 2 . These operators are well-defined for w G ^5,^.0 
so long as the system Eq. (2.1) is a linear symmetric hyperbolic Fuchsian system for 
parameters 5, fi, r, q and go as in Definition 2.5. We now define w to be a weak solu- 
tion of Eq. (2.1) with vanishing leading term provided it satisfies, for all test functions 

0, 

(r[w}^):=(C[w}-T[w},4>)=0. (2.25) 

Here, we note the discussion of distributional time derivatives in Section A of the ap- 
pendix. 

Before proceeding to show that weak solutions exist, we establish the following useful 
technical result. 

Lemma 2.9. Suppose that Eq. (2.1) satisfies the conditions to be a linear symmetric 
hyperbolic Fuchsian system for a fixed set of parameters 5, fi, (, 5, r, q, go as per 
Definition 2.5, and is block diagonal with respect to \i. Then for every test function <j), 
the maps {£[■], (ft) and {J~ [•],(/>) are bounded linear functionals on X$^£. 

Proof. To prove this lemma it is sufficient to show that each term in 4>) is bounded 

by C\ \w\ \s,^,o, for some positive constant C and for every w G ^5,/i,o- We demonstrate 
this for the first term, J Q (TZ[fj]Siw, D(j>) L2 dt. Using Holder's inequality, the spatial 
continuity 3 of Si and the block-diagonal property, we find that 



I 



8 

{K[fi]Siw,D<l>) L 2dt 



Jo 

<5 sup I |7i[/i]5iu;||i3| p^Hxa < CIMk^o- 
te(o,<5] 



x The operator L is the adjoint form of the linear principal part operator L (see Eq. (2.13)). 

2 Note that we do not need to set <^>(0, x) = in order to avoid boundary terms at t = 0; the decay of 



w(t, ■) and of /o(£, •) as t approaches zero removes such potential boundary terms 

3 This fo 
embedding 



3 This follows from the definition of a linear symmetric hyperbolic Fuchsian system, and from Sobolev 
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The constant C, which is used to estimate both the contributions from S\ and from <f>, 
is uniform in the sense defined above. Other terms in (£[u>], (j>) follow similarly, and the 
same arguments hold for the {J-[-},4>) operator. □ 

We now determine that, for a given linear symmetric hyperbolic Fuchsian system 
with certain conditions holding, the singular initial value problem with zero leading 
term has a weak solution. In the proof of this result, we show that these solutions can 
be obtained as a limit of approximate solutions of the Cauchy problem, as described 
above. 

Proposition 2.10 (Existence of weak solutions of the linear singular initial value prob- 
lem with smooth coefficients). Suppose that Eq. (2.1) satisfies the same conditions as 
stated in Lemma 2. 7 and hence is a linear symmetric hyperbolic Fuchsian system ( with 
smooth coefficients) for 5, fx, (, r, q, and qo as per Definition 2.5 with q = and qo = 2. 
Then there exist weak solutions w : (0, 6] x T 1 — > M. d to the singular initial value problem 
(with vanishing leading term) which are elements of X$ jfij o. 

Observe that Proposition 2.10 is the most general existence result which we obtain for 
linear equations, in the sense that only minimal control of the behavior of the coefficients 
of the equation is required (i.e., q = 0, qo = 2 as in Lemma 2.7). We discuss higher 
regularity of the solutions under stronger regularity assumptions in Proposition 2.12 
below. We also note that while Proposition 2.10 provides sufficient conditions for the 
existence of solutions, it tells us nothing regarding uniqueness. To obtain uniqueness, 
we need to impose stronger assumptions on the coefficients of the PDE system; see 
Proposition 2.14 below. 

Proof. As described above, we choose a sequence (t n ) converging to zero, and the corre- 
sponding sequence of approximate solutions (v n ) € C°((0, 5], H qo (T 1 )) n Xg^ Q. We 
seek to show that the sequence (v n ) forms a Cauchy sequence in Xg^ t Q. Defining 
v m - v n , we readily see that 

r 0, t E (0,t m ], 

£,mn(t)%) — \ V mi t £ {trnitri\i (2.26) 

[ v m -v n , t 6 (t n ,S\. 

From the energy estimate for the Cauchy problem Lemma 2.7 on each subinterval, we 
then derive 

r = o, te(o,g, 

HW(V)fmn(V)IUM < O + Cj^s^WRMfoWvds, t£(t m ,t n ], (2.27) 

I < \\TZ[fi](tZ-)v m (t n ,-)\\ L2 , te (t n ,S\, 

where in the last inequality we have used the energy/norm equivalence Eq. (2.15) above, 
and we have also used the fact that the (linear) PDE system for v m — v n has a vanishing 
source term /q- Recalling the definition of the norm || • | \s,(i,q, noting the monotonicity 
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of f* s 1 1 |7^[/i]/o | \\ ds, and noting the equality £, mn (t n , 
have 

1 1 U 1 1 8,11,0 = sup \\K[n](t,-)£ mn (t,-)\\ L 2 <C 

te(o,8] 



•) = v m (t n , •) for t = t n , we now 



s 1 ||^M/o||x,2(is. 



To complete the argument that we have a Cauchy sequence, it is useful to introduce 



G(t) := / s- 1 \\TZ[ij}fo\\ L 2(s)d S , 



(2.28) 



which is well-defined so long as /o € Xs tV ,o for v > fx. Choosing e > as a lower bound 
for the gap between v and fx among all components, we see that there must exist a 
constant C such that G(t) < Ci e ; thence, we have 



I |£mn| |<5,/i,0 

< C\G(t n ) - G(t 



(2.29) 



from which it easily follows that (v n ) is a Cauchy sequence in the Banach space Xg^Q. 

Since it has been established (in Lemma 2.9) that V = C — J- is a continuous operator 
on Xs,^,o, to show that the limit of the Cauchy sequence (v n ) is a weak solution of 
the system of interest, it is sufficient to show that the limit of the sequence of reals 
((^bn]) <j>)) is zero for all test functions (f>. Choosing any v n in our sequence, we know 
from its definition that v n vanishes on (0, t n \ and is a solution to the equation {V[-], 4>) = 
on [t n , 8]. Recalling the definition of V, we calculate on this latter interval, for any test 
function 6. 



\{V[Vn]A)\ 



{K[lAhi<t>)i?(Tt) 



dt 



Straightforward calculation then shows that 



o> Wz,2( T i) 



dt < 



< 



(^•M/0»L2 (T i) \dt 



dx\n[fi]f \ 2 



1/2 



dx\n[ii\f \ 2 



1/2 



dx\^\ 2 ) 
dx\4>\ 2 

T 1 



dt 

1/2 



L2 / r l WR[n}Ut)\\ L 2dt<CG{t n ) 



dt 



< sup 

te(o,8] 

from which it follows (from the properties of G(t)), that we have a weak solution. □ 

Based on this existence result for weak solutions, we would like to define a map 
which, for a fixed choice of Si, S2, N and F\, maps any smooth function fa S Xs, u ,o t° a 
weak solution w G Xg v $ of (V[v n ], 4>) = 0. Then as a next step, we would like to extend 
this map to all /o of Xs, u ,o, and thereby show that weak solutions exist for all /o £ Xg, v ,o , 
and not just for those /o which are smooth. While the lack of a uniqueness result for 
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weak solutions is an impediment to defining the desired map, we can get around this by 
provisionally defining an operator of this sort which maps a smooth fo to the solution 
of the weak solution which is obtained as the limit of the sequence (v n ) (as discussed in 
the proof of Proposition 2.10). We do this now, noting that the definition makes sense 
only after we have established that we get the same limit for any choice of the sequence 
of times t n ). We then establish an estimate for this operator, and use this estimate to 
extend the operator to all of Xg^Q. 

Proposition 2.11. Presuming the hypotheses listed in Proposition 2.10, there exists an 
operator H : X§ v q —¥ Xguo which maps a smooth source function fo to the weak solution 
w of [Vim], (ft) = which is obtained as the limit of the sequence of approximate solutions 
(v n ) corresponding to a choice of a monotonic sequence of times (t n ) converging to zero. 
This operator is well-defined (independent of the choice of the sequence (t n ) ) and satisfies 
the estimate 

||H[/o]||^,o<^||/o||^,o, (2.30) 

for all smooth fo £ X$ v q. The positive constants C and p are uniform. 

The operator extends to all (not necessarily smooth) fo £ Xg^o, with the estimate 
(2.30) holding for all such fo with the same constants. Indeed, this extended operator H 
maps all fo £ Xg v $ to weak solutions of Eq. (2.1). 

The last paragraph in this proposition generalizes the existence result in Proposi- 
tion 2.10 to all, not necessarily smooth, source terms fo £ Xg iU; o- We note, however, 
that otherwise the system is still assumed to have smooth coefficients in the sense defined 
above. 

Proof. We presume initially (as part of the hypothesis of smooth coefficients) that fo is 
smooth; i.e., f £ C°°((0,<5] x T 1 ) n X s , u ,o- To show that H is a well-defined map from 
C°°((0, 5] x T 1 ) n Xg Vj Q to Xs :fJ ,fii independent of the choice of time sequence, we choose 
a pair of such sequences (t^) and (i^J with their corresponding sequences (v^) and (v^) 
of approximate solutions, and from the union of the two time sequences we construct 
a third time sequence (ti). As is the case for (v\) and (v^), the combined sequence of 
approximate solutions (vi) must be a Cauchy sequence, so 1 — f^JI<5,/i,o must vanish 
in the limit n, m — > oo. Then labeling w 1 as the limit of the first sequence and w 2 as 
the limit of the second, we calculate 

Ww 1 ~ ^ 2 |k^o < Ik 1 - vlWs^o + \\vfn - w 2 \\ s ,n,o + \\vl - ?4ll<W)- 

It easily follows that w 1 and w 2 are equal in X$^$- 

To prove the estimate for H (restricted to smooth fo), we let (v n ) be a sequence of 
approximate solutions with limit w = H(/o), and then based on Eq. (2.29) we determine 
that \\w - < CG(h) < CG(5). It then follows that 

Ikll^o < IN \\s,fi,o + CG(S). 

1 Here, we set S to be the smallest bound among the two sequences. 
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If we now apply the energy estimates to show that H^iH^o < CG(5), we deduce that 

IMk/i,o < CG(8), 

for some adapted constant C. To relate G(S) to the source term, we check that 

s- p \\n^}f \\ L 2 < ||/oll«w) 

for some p > so long as p < v. It then follows from multiplying both sides by s _1 and 
integrating over that 

G(S) < i^H/oll^o- 
P 

The estimate Eq. (2.30) is then a consequence. 

To extend the domain of H from C°°((0, 6] x T 1 ) n X 5)i , )0 to X«j )1/0 , we note that this 
first space is dense in the second by definition. Hence, for any /o G Xg V Q, we can find 
a sequence of functions /oj G C°°((0, 6] x T 1 ) n Xg ;U ,o which converges to /o- It follows 
as a consequence of the estimate Eq. (2.30) that there is a unique continuous extension 
of H to the full space Xg v o- The extended operator, which we refer to with the same 
symbol H, is continuous and satisfies the same estimate. 

The continuity of the extended operator H and the continuity of {V [w] , <j>) = with 
respect to w easily implies that EI maps any fo G X§ v q, even those which are not smooth, 
to weak solutions. □ 

To proceed from weak solutions to strong solutions of the singular initial value prob- 
lem for these linear systems (while still keeping the smoothness assumption for the 
coefficients Si t i, S2,l> -^l and F\), we need to determine the regularity of these weak 
solutions. We do this in the following proposition, and thereby prove the existence of 
strong solutions. 

Proposition 2.12 (Regularity of solutions for smooth coefficients). Suppose that all of 
the conditions of Proposition 2.10 hold, with the exception that q > 1 and qo = q + 2. 
Then, weak solutions w of the singular initial value problem (whose existence has been 
checked in Proposition 2.10) are differentiable in time 1 and hence are strong solutions of 
Eq. (2.1), with w G and Dw G X$ ^q—l- As well, the solution operator H defined 

in Proposition 2.10 maps X$ vq to X§ „ g , and satisfies 

WmMhw < 6 p C\\f \\ St „ >q , (2.31) 

for all (not necessarily smooth) fo G X$ uq . The constants C > and p > are uniform 
in the sense of Definition 2.6 (but may depend in particular on q). 

Observe (without pursuing the details here) that an estimate similar to Eq. (2.31) 
can also be proven for the time derivative of the solution. Additional regularity assump- 
tions on the time derivatives of the coefficients of the equation also allow one to prove 
corresponding statements regarding higher order time derivatives D k w for k' > 2. 

in the sense of Section A in the appendix 
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Proof. Using w to denote the solution to the singular initial value problem whose ex- 
istence is established in Proposition 2.10 (as an element of Xg „o), and using (v n ) to 
denote the sequence of approximate solutions which converges to w, we first note that 
it follows from their definitions 1 that the v n are contained in C°((0, 6], H q (T 1 )) - even 
in Xs tflt q - under the hypothesis of Proposition 2.12. Hence, in the same way as we have 
used the energy estimates in Lemma 2.7 to show that (v n ) is a Cauchy sequence in Xg „( 
for the proof of Proposition 2.10, we can now use the energy estimates of Lemma 2.8 to 
show that (v n ) is a Cauchy sequence in Xg jfJ/ - e ,q for an arbitrarily small constant e > 0. 
We then show that the limit of this Cauchy sequence equals w above. The solution w is 
hence in Xg^_ 6iq and we get the estimate 

\\nfo]h n <s p c\\f \u,v, q . 

Now, if the equation is of linear symmetric hyperbolic Fuchsian form for a choice of 
/i, as we have assumed so far, it is also of linear symmetric hyperbolic Fuchsian form for 
a choice of /I := /i + e if e > is sufficiently small in comparison to v — \l. Moreover, the 
assumption that the system is block diagonal with respect to ji also implies that this 
is the case with respect to /2; the same is true for the condition involving the energy 
dissipation matrix. Hence, we can apply the argument in the previous paragraph based 
on this choice of ju. This leads to the conclusion that, in fact, the solution w is in 
x Sfi-e, q = x &,n,q ( as opposed to Xfi^-e^ above) and 

P[/o]lk M , g <^C||/ |k^, 

possibly after a slight change of the constants C and p. 

Next we show that the solution w is differentiable in time. We define 

v n := Sf 1 (-S 2 td x v n - Nv n + f + F x v n ) . 

We know that v n G Xg^ q _i and v n (t) = Dv n (t) for all t G [Si, 5] for any 5j G (0, 5) and 
for n sufficiently large; we cannot choose 5i = here since the time derivative of v n is in 
general not defined at t = 0. Moreover, we find from the definition and the convergence 
of the sequence v n that 

\\ v n ~ v m\\&,ijL,q— 1 £: C||t>n — V m \ \ g^q 

for a uniform constant C > 0. Hence there exists v G Xg„ q _x such that v n — > v. The 
estimate also holds if we restrict the time interval to [5i , 5] as above and hence we find 
that Dv n (t) = v n (t) — > v(t) uniformly at every t G [5j, 5]. It is then a standard result 
that w is differentiable in t at every t G [Si, 6] and further that v(t) = Dw(t). Since 8i 
can be chosen arbitrarily small, it follows that for all t G (0,5], w is differentiable in t 
and v = Dw. Consequently, we find that Dw = v£ Xg^ q -\. 

To argue that w is a strong solution, we start from the fact that w is a solution of 
the weak equation whose integral representation can be integrated by parts in both time 

J By definition, v n (t, •) = for t 6 (0, t n ], and for t £ [t n , 5] it is equal to the unique solution of the 
Cauchy problem with zero data at t n . 
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(using Eq. (A. 4) in the appendix) and space. We may then choose a suitable sequence 
of test functions, for example those which are used as mollifiers in Lemma A.l in the 
appendix, so that the resulting system converges pointwise almost everywhere to one of 
the components of Eq. (2.1) evaluated at one point (t, x). Doing this for every component 
and for every point (t, x) G (0, 6] x T , we determine that w is actually a solution of the 
strong equation almost everywhere. □ 

To this point, we have assumed throughout our analysis that the matrices Si i, S2,i, 
Ni and F\ are smooth; i.e., we have not thus far allowed these matrices to be general 
elements of the spaces Bg ^ q r from Definition 2.5. If we wish to use our current (linear) 
results as a tool for proving that there are (unique) solutions to the singular initial value 
problem for the (nonlinear) quasilinear system, we need to generalize these linear results 
to include the possibility that the matrices listed above are not smooth (since, in the 
quasilinear case, these matrices are functions of the solutions, which may not a priori be 
smooth). 

Before carrying out this generalization of the existence (and uniqueness) results for 
the linear singular initial value problem, we note that we can at this stage assume that the 
term Fx vanishes. This simplification does not constitute an essential loss of generality 
since in our work below, we replace the linear source term function /o by a general 
quasilinear expression shortly; the resulting expression then incorporates the effects of 
the term F\. We recall that the Fi-term plays a convenient role in our verification that 
the higher order energy estimates of Lemma 2.8 hold. Such a term is generated in the 
linear equation for d x v, which we obtain by taking a spatial derivative of the equation 
for v. 

Proposition 2.13 (Existence for the linear singular initial value problem for non-smooth 
coefficients). Suppose that Eq. (2.1) is a linear symmetric hyperbolic Fuchsian system for 
5, p,, £, q, qo and r as in Definition 2.5 (with F\ = 0) for qo = q + 2 and q > 2 (not 
necessarily with smooth coefficients), and suppose that it is block diagonal with respect 
to p. Suppose that the energy dissipation matrix Eq. (2.6) is positive definite. Then 
there exists a solution w : (0, 5] x T 1 — > M. d to the singular initial value problem with 
zero leading order term such that w G Xg^q and Dw G Xs^^-x- The solution operator 
H : fo I—)- w maps Xg^q into Xs^ u ,q, and satisfies 

P[/o]|| 5 , M <5 p C||/o|| 5 ,, i(? 

for some positive uniform constants C and p. 

Observe that this result also holds in the case that qo and q are both infinite: If the 
conditions of this proposition are satisfied for all integers qo = q + 2 and q > 2, then 
w G X s . /i,oo and Dw G X§^ >00 . However, the ^-parametrized sequence of constants C 
and p occurring in the estimate of the solution operator may in general be unbounded 
as q — > oo. 

Proof. The basic idea is to approximate the non-smooth coefficients Si x, S2,i and N\ by 
a sequence of smooth ones and then to apply Proposition 2.12 to obtain a sequence of 
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approximate solutions. The main work is then to prove that this sequence converges to 
the solution of the system with non-smooth coefficients in an appropriate sense. 
Step 1: A sequence of approximate solutions. We presume that a linear symmetric hyper- 
bolic Fuchsian system with parameters 5, r, q, qo, \i and £, and with coefficient matrices 
5i,i, £2,1 and Ni in the perturbation space Bs^,q,r has been specified. We assume that 
this system is block diagonal with respect to \x and that the energy dissipation matrix is 
positive definite. According to the definition of the space B$£^ r , there exist sequences 
0^2,i,h']) and (^l.b'l) °f smooth elements in Bg^ iqr which converge to 51,1,52,1 
and N\ (in a way described below). We thus obtain a sequence of linear principal part 
operators (with smooth coefficients) 

L b] [w] := (5i, + S hm )Dw + K[l - /3](5 2 , + S 2>m )td x w + (N + N m )w, (2.32) 

and hence a sequence of systems of the form Lmfto] = /o- For each j, this is a linear 
symmetric hyperbolic Fuchsian system for 5, fi, q, qo, £ and r with smooth coefficients. 
It is clear that the sequences can be chosen so that, for each j, the block diagonal 
condition with respect to [i is satisfied, and the energy dissipation matrix is positive 
definite for each equation. Clearly as well, each 5i,i ui is differentiable in time and 
DS\ \ ui is bounded (in a sense which we make more precise below). Hence, for each 
j, Proposition 2.12 implies the existence of a solution operator Hm, and therefore a 
sequence wyj G -^S,fj,,g defined by wyj := H^jf/o]. 

Step 2: Uniformity of the sequence of coefficients. To study the convergence properties of 
these approximate solutions wui in more detail we make a special choice of the sequences 
(5iim), OS2 im) and (JViu) as in Lemma A.l in the appendix (where we replace the 
two indices i and j by just one index j). The advantage of this choice is that we will 
be able to argue that \\1Z[(]DSi^^j](t, •)\\[ J 2 is uniformly bounded in j and t under the 
hypotheses of Proposition 2.13, which will be important for the following argument. 
The slight disadvantage, however, as stated in Condition (iv) of Lemma A.l, is that the 
convergence is guaranteed only with respect to a norm || • |L-r for any exponent matrix 

C smaller than £; fortunately we will see below that this is not significant. Let us choose 
such an exponent matrix £ with strictly positive entries. Let us moreover suppose for the 
moment that a uniform bound for ||7?.[^]Z)5i i i u-i (t, -)\\l 2 has been found (which we show 
shortly) . By setting £ = £ and choosing some uniform value of the (function space ball) 
radius s in Condition (ii) of Proposition 2.12, we are allowed to apply Proposition 2.12 
in such a way that each of the approximate equations Luj [w] = fo is a perturbation of 
one common equation in the perturbation space B s ^ r . A particular consequence is 
then that we obtain an estimate for the operators Hm of the form Eq. (2.31) with C 
independent of j. 

To establish this uniform bound of ||72{£].D5i ) i ) u](t, -)IIl 2 > we use Eqs. (A.l) and (A. 2) 
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from Lemma A.l in the appendix to obtain 



K[C]DS ltm =Ki(]DKi-(]Ki(]S ltUj] 



+ 



(W[C]5i,i)(fl,y)— 



{-i)—t4> 

a- 

3 



a ; 



dy ds, 



where 5i i is the extension introduced in Lemma A.l. We wish to estimate this expres- 
sion in the L 2 -norm. The first term can be estimated in the L 2 -norm by C(^)||S , 11 jj] 

with a constant determined by £. The second term is treated as follows 1 
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where we have used the definition of the extension Si i in the second line. The properties 
of the kernel cf> imply that the term with the L^-norm is unity (independently of x). Since 
q > 1, we can use Sobolev embedding to estimate the term with the L^°-norm by the 
Hy-noim and hence the sup s -term by the norm || • II5 f 9 - As a consequence, all quantities 

are independent of x, and therefore we find that the second term of lZ[QDSi t iu-\ above 
can in total be estimated in the L 2 -norm as less than or equal to 
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for some constant C > which is, in particular, independent of t and j. Hence, we only 
need to estimate the last integral: 
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This takes the form x ■ <f>(x) for all This has the property that it vanishes for all 

\x\ > 1 and is bounded for all |x| < 1. Summarizing, we have thus confirmed that 
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1 We write e.g. L^iT 1 ) to fix the integration variable corresponding to the norm. 
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for every j, and for a constant C, which, in particular, does not depend on j. 

Step 3: Convergence of the sequence {wy\ ) in Xs lfJ ,,q—i- We seek to show that the sequence 
(wyi) converges 1 in Xg^ q _i. We do this by showing that (wyj) is a Cauchy sequence: 
Setting := — uuyj, we derive the equation 

L [i\C[ij] = ~( L [i\ ~ L \j])w[j], (2.33) 

where we interpret the right hand side as a source term for this linear equation for £^1. 
One readily checks that Eq. (2.33) is a linear symmetric hyperbolic system for the same 
parameters as above, but with differentiability index q — 1 (since the source term incor- 
porates one spatial derivative); hence, so long as q > 2, we may apply Proposition 2.12 
also to this equation. We thus obtain, as a consequence of Eq. (2.31), 

IK] - w\j]\kn, q -i < c - £[j]Mj]|| v+ j mini(; _i , 

where the (scalar) constant Cmin is the minimal value of all components of £ at all x € T 1 
(note that Cmin is positive), and where the generic constant C represents the constant 
in Eq. (2.31) (for v replaced by \i + Cmin)- It is crucial here that the constant C does 
not depend on the index i; this is a consequence of the uniformity of the constant in 
Eq. (2.31). If we now expand out the definition of Ly\, using the block diagonality 
conditions, the Sobolev embedding (for spatial dimension one and q > 2) and the Moser 
inequality (stated for instance in Proposition 3.7 in Chapter 13 of [40]), we find that 

IK] - wyjik^g-i < c(||<Si,i,m - s i,mh,iq-i + W s w,i<i ~ s 2,mh,iq-i 

+ \\N li{ n - N m \\ 5 ^ q _^j \\fo\U,u,q- 

It immediately follows that (wy^) is a Cauchy sequence in Xs^ q -i, and hence the se- 
quence has a limit w £ Xg 

Step 4-' The limit w is a solution of the original equation. Standard arguments of the 
sort carried out in the proof of Proposition 2.10 show that w is a weak solution of the 
system Lw = /o with non-smooth coefficients. Since each wyi is differentiable in time 
and a strong solution of the equation, we can solve each equation for Dwyj . Relying on 
an argument similar to that used in the proof of Proposition 2.12, we see that it follows 
that w is differentiable in time, with Dw 6 X$^^ q _2, and therefore the limit w is actually 
a strong solution of the equation. 

Step 5: The limit w is in Xg^ iq . We now show that w is in fact in X$^ q (and not 
only in Xs^^-i) and consequently Dw S Xs^^-i (and not only in Xs^ )Q -2)- It follows 
from its definition that the sequence (wyj) is contained in Xg^g, and furthermore, as a 
consequence of the operator estimate Eq. (2.31), we have that 

||«>[j]lk M) g < C 1 1 /o|| S,v,q- (2-34) 
x The reason for considering q — 1 and not q is made clear below. 
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Here, the constant C is independent of the index j as a result of our discussion of 
uniformity above. We thus find that the sequence is uniformly bounded in Xg a,q- If we 
now fix a time to £ (0,5), then the sequence (w^(to)) is bounded in the Hilbert space 
-ff^T 1 ). Since the sequence (w^(to)) is also convergent in the Hilbert space H q ~ l {T l ), 
we can apply Corollary C.3 from the Appendix. We hence find that the limit w(to) 
is in H q (T l ). We consider the function w, which is a strong solution in Xg t u q -i with 
w(to) £ H q (T 1 ), to be a strong solution of the Cauchy problem of the linear symmetric 
hyperbolic equation with the "initial data" w(to) contained in H q (T l ). Given that the 
coefficients of the system are contained in C°((0, 5], H q (T )), the standard theory of 
linear symmetric hyperbolic equations (see [40]) implies that w, and therefore lZ[[i]w 
(recall that K[n] is in C°°((0,<5] x T 1 )) is a continuous map {0,5] -> ^(T 1 ). In fact, 
this latter map is bounded as a consequence of taking the limit j — > oo of Eq. (2.34); i.e., 
w £ Xs^ q (see the appendix). Replacing \i by \x + ?in all of the previous steps, we see 
that the same arguments go through as long as e > is sufficiently small. We therefore 
find that TZ[fi + e\w is a bounded continuous map (0, 6] — > H q (T 1 ). It then follows from 
Lemma A. 2 in the appendix that w £ Xgnq- 

Step 6: Properties of the solution operator H. We have thus extended the solution 
operator M to the case of non-smooth coefficients by the above limit procedure. It can 
thus be checked that the estimate Eq. (2.31) still holds with uniform constants. □ 

Having obtained a comprehensive existence result for the singular initial value prob- 
lem of linear symmetric hyperbolic Fuchsian equations, we now show uniqueness of these 
solutions. 

Proposition 2.14 (Uniqueness for the linear singular initial value problem). Suppose 
that all of the conditions of Proposition 2.13 hold for a chosen singular initial value 
problem (with zero leading-order term). The solution for this problem is unique in Xg^q- 

Proof. We consider w and w to be a pair of (generally different) solutions to the same 
singular initial value problem, and we define oo := w — w to be the difference between the 
two. It follows that w is a solution of the same equation with vanishing source-term /o, 
with ui(t) being an element of H 2 (T l ) for every time t G (0, 5]. Choosing any to £ (0, 5], 
we can also consider w|(t 0)( 5] to be the unique solution of the Cauchy initial value problem 
(for the same linear PDE system) with initial data u(to). Since the solution uj together 
with the coefficients have i7 2 -regularity and since S\ is guaranteed to be positive definite 
on the whole time interval, we may apply the energy estimate Eq. (2.17). We obtain 
(replacing \x by \x — e, which is allowed for any e > 0) 

\\U\p - e]u\\ L 2(t) < C\\K[ii - e]uj\\ L 2(t ), (2.35) 

for all t £ (to>^]i with the constant C independent of t. Observe that it follows from the 
definition of H[fi] that we can rewrite the right hand side of Eq. (2.35) as 

C\\K[fi - e}u\\ L 2(t ) = ^CH^HI^o), (2.36) 
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so that Eq. (2.35) takes the form 



\\n[fi - e]oj\\ L2 (t) < tlC\\K[ti}u\\ L2 (to). (2.37) 

We now take the limits of Eq. (2.37) as to ~~ * 0, noting that the left hand side of the 
equation and the constants are unchanged by taking this limit. Since lZ[fj\u is a bounded 
map from (0, 5] to L 2 (T 1 ), the limit as to of the right hand side of Eq. (2.35) vanishes. 
It thus follows that for all t 6 (0, 5], 

\\TZ[fx-e]co\\ L2 (t) = 0. 

Then since 1Z[fi — e](t,x) is bounded positive at any fixed t on (0,5], we deduce that 
u(t,x) = at all t; uniqueness follows. □ 

2.3.3 The nonlinear theory 

The results obtained in Section 2.3.2 pertain exclusively to linear systems. In this section, 
we use those results together with a fixed point iteration procedure to prove Theorem 
2.4, which establishes existence and uniqueness of solutions to the singular initial value 
problem for the nonlinear system Eq. (2.1) with a (no longer necessarily vanishing) 
leading-order term uq. 

To start, it is useful to rewrite Eq. (2.1) in a convenient form. Recalling the definition 
(see Eq. (2.5)) of the operator L(u)[v] := Si(u)Dv + S2(u)td x v + N(u)v, we may write 
Eq. (2.1) in the form L(u)[u] = f(u). Despite the nonlinear nature of L, this operator 
is linear in the sense that L(u)[v\ + V2] = L(u)[v\] + L(u)\v2[. Hence if we let no 
denote a chosen leading-order term (satisfying the hypotheses of Theorem 2.4), if we set 
u = uq + w, and if we recall the definition (see Eq. (2.7)) 

&(uo)[w] := F(it )H ~ L(u + w)[u ], 

then Eq. (2.1) takes the form 

L(u + w) [w] = 3"(n )H . (2.38) 

The linear analysis discussed in Section 2.3.2 does not apply to the nonlinear equation 
Eq. (2.38) directly. Nevertheless, if we linearize this equation by fixing w £ Bs tfJl ,q tS (for 
some s > 0) and write 

L{u + w)[w] = ?(u )[w\, (2.39) 

then (presuming the hypothesis of Theorem 2.4) the techniques of Section 2.3.2 are 
applicable. In applying these techniques, we assume that 5 has been chosen sufficiently 
small so that Condition (ii) of Definition 2.5 is satisfied; this implies no loss of generality 
since, as we see below, the argument leading to the proof of Theorem 2.4 requires further 
shrinkage of the time interval. It is important to note that, for every w G B§„ qs , the 
hypothesis of Theorem 2.4 implies the existence of common quantities £ and r (as in 
Definition 2.2) so that Si t i(uo + w), S2,i(uq-\-w) and Ni(uq + w) are all contained Bg^ iqr . 
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Moreover, if to € Bg/ jfiqjS for any 6' < 5, then the same statement regarding Sx,ij ^2,1 
and Ni holds for the same common £ and r, but with 5 replaced by 5' . 

Replacing the right-hand side of this equation by a fixed function 4> G Xs tU>q , we 
readily check that the linear system 

L(u + w)[w}(t,x) = 4>(t,x), (2.40) 

is of linear symmetric hyperbolic Fuchsian form (Definition 2.5) for q and qo = q + 2 and 
for a sufficiently large constant r. Hence, it follows from Propositions 2.13 and 2.14 that 
the system has a unique solution w G Xg^ q (we only require q > 2 at this stage of the 
proof), and we can define the corresponding solution operator W(uq + w) which maps 
the source term cj> to the solution w = M(uq + w)[4>]. The case cj> = 3^(uq)[vj\ corresponds 
to Eq. (2.39); thus we compose M(uo + w) with 3~(«o) to define the operator G(uq) as 
follows: 

G{u )[w] := M{u + w)[?(u )[w\]. 

Hence, w = &(uq)[w\ G Xg ifJLi q is the unique solution of the singular initial value problem 
of Eq. (2.39). In terms of G(uq), we see that w is a solution of the singular initial value 
problem for the nonlinear equation (2.1) with leading-order term uq if and only if it 
satisfies w = G(uq)[w]; i.e., if and only if w is a fixed point of G(uq). 

The operator G(uq) is the key to the following fixed point iteration argument. We 
define the sequence of functions (w n) by setting u>o = 0, and defining w^ + i = G(uq)[wn] 
for N G N. To control this sequence, we need uniform bounds; i.e., we wish to show that 
each element of the sequence is contained in Bg tfl ,q,s- Suppose that this is true for 
Wo, . . . , u>n- It follows from the hypothesis of Theorem 2.4 (given that wq = 0) that 

muo)[w N ]\\ s , u>q < C\\w N \\ Sjm + \\?(u )[Q]\\ 5)U! q. 

The constant C > does not depend on N. Using the definition of iojv+1, together with 
Eq. (2.31), we have that 

\\wN+i\\8 4 i, q < 5 p C||J(« )Kv]lk!,,g' 

where C and p > are constants which also do not depend on N. Combining, we obtain 

\\w N+1 \\ s ^ q < 6PC(C\\w n \\sm+ Wuo)[0]\\s, v , q ) ■ 

We recall that the uniformity of the constants implies that the same estimate holds with 
the same constants if we choose to formulate the same singular initial value problem in 
terms of a constant 5 G (0, 5) instead of S itself. Since we have supposed that wn is 
contained in Bg^q^ ~ that is, jv||<5,/Lt,q < s - we can find such a sufficiently small 5 so 
that 

l p CC < 1/2 and ^CH^UoMO] \\s,u,q < s/2, 

while preserving the bound ||i«jv||a«g < s. This can be done, since ||3 : '('Uo)[0]||| ^ g < 
||3 ! '('Wo)[0] ll^j/,^- For this diminished choice 5, we thus determine that ||iojv+i|||„ q < s. 



33 



Since the above estimates do not depend on the index N, it follows that the whole 
sequence is bounded, and we have (wn) C B$ n q s - 

We now consider an arbitrary pair of functions w,v £ Bg n qs , and we calculate the 
following estimate for the norm of the difference of the operator G(uq) acting on each of 
these: 

||G(«o)[tu] - G(^o)Mll5 )At , g _i < P(u + w)[?(u )[w}} - M(u + ^[^oMIIj^-i 

+ ||H(«o + iu)[2H)M] -H(«o + «)[3-(uo)H]|| ?lM ,,-i. 

Note that, for reasons discussed below, we work with II • lie „ „_■. rather than II • II? ,, „. It 

ii nu,/.t,t/ _l ii no,/.t,y 

follows from the hypothesis of Theorem 2.4 and from Eq. (2.31) that the first term on 
the right hand side of this estimate satisfies the inequality 

||H(«o + w)[?(u )[w}] - H(u + ^)^(«o)H]|| 5 >, 3 _i <C\\w- v^^, 

with the Lipschitz constant C smaller than unity so long as we allow a further decrease 
in 5; the argument for this is the same as for the semilinear case [13, 12]. This controls 
this first term. 

To estimate the second term on the right hand side, we set 

w A ■= H(n + w)[2(uo)[v]], w B := H(u + v)[$(uq)[v]]. 

It follows from the definition of M(uq + w) that 

L(u + w)[w A ] = 3"(n )M, L(u + v)[w B ] = &(uo)[v]. 

Therefore, setting L(u$ + w)[wa] and L(uq + v)[wb] equal, and using the linear property 
of the operator L(uq + w) noted above, we derive 

L(uq + w)[wa - w B ] = -(L(uq + w) - L(u + v))[wb]- ( 2 -41) 

The right hand side of this equation is similar to a term which appears in Eq. (2.33); 
thus we can treat it using similar techniques to those used in the proof of Proposition 
2.13. In doing this, we rely on the hypothesis for Theorem 2.4, and we use the condition 
q > 3 in order to guarantee that the source term of Eq. (2.41) has at least two spatial 
derivatives. We thus obtain 

||M(«o + w)[F(u )[v\] - M(n + ^MMlll^-i 
< 5Pd(\\S 1A (w) - S 1 , l (v)\\- S ^ 1 + \\S 2tl (w) - ^,i(«)||j )C , ff -i 
+ \\N 1 (w)-N 1 (v)\\ Ujq _ 1 ) 

for a constant C which may depend on s, but not on the particular choice of 5. 

Again using the hypothesis of Theorem 2.4, we see that for a choice of a possibly even 
smaller 6, which we now label 6, we can control this second term from the right hand side 
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of the estimate for ||G(iio)[to] — G(?xo)[f]|| j q _ l via a term of the form C||u> — v\\g t for 
C € (0, 1). We thus determine that indeed the operator G(tto) is a contraction mapping 
on -Bj^ 9 _ ls (for sufficiently small 5). It follows from standard arguments that the 
sequence (wn) has a unique limit w, contained in Bg _^ s , which is a fixed point for 
G(«o) and hence is a weak solution. 

The sequence (wn) C Bg ^ g is bounded in Xg , but to this stage is known only to 
converge in Xg ^ to w. This situation is similar to that encountered in the proof of 
Proposition 2.13. A similar argument involving Corollary C.3 and the standard Cauchy 
problem of hyperbolic equations implies that w is indeed an element of Xg . 

To show that w is the remainder of a strong solution of the singular initial value 
problem, it remains for us to check that w is differentiable in time. The definition of the 
sequence (wn) shows that for each integer iV, Dwn exists and is contained in Xg v 
Furthermore, this sequence converges in Xg _j by a similar argument as in the proof of 
Proposition 2.12 using the Condition (iii), and the positivity of £. Since this convergence 
is uniform in time, it follows that w is differentiable at each t and that Dw(t) is the limit 
of (DwN(t)) at each t. It follows from this limiting procedure that u is indeed a strong 
solution to the singular initial value problem; similar arguments have been used before 
also in the proof of Proposition 2.13. 

The uniqueness of this solution w follows from the uniqueness of the fixed point for 
the contraction mapping G(ito). 

In order to complete the proof of Theorem 2.4, we must consider the case q = oo. 
We do this inductively in q. It is important to notice here that all of the constants in 
the previous estimates may depend on q and hence we may have to adapt the choice of 
5 in each induction step. It is thus possible that the sequence of these constants (8 q ) 
tends to zero as q — > oo. To show that this possibility is avoided we use the result that 
any solution to the Cauchy problem for symmetric hyperbolic systems with a bounded 
first spatial derivative can be extended to a common time interval. Let us fix any 
q > 3. Theorem 2.4 (with finite q) shows that there exists a solution w G Xj^^ for 

some 5 G (0,(5]. Let i* G (0,5], and consider the regular Cauchy problem with data 
w(t*,x). Since <7 ^2, the Sobolev inequalities guarantee that the first spatial derivative 
is bounded on Ji*, 5]^and we may apply Proposition 1.5 in Chapter 16 of [40] to show that 
there exists a 62 > 5 such that the solution may be extended as a i/ 9 -solution to (0, 62]- 
The same argument applied to any other value of q > 3 implies that the solution can be 
extended as i^-solutions to the same time interval (0, 62]- For q = 00, we therefore find 
a unique solution w on the same time interval in Xj (and hence Dw G X-? ). 

2.4 Existence and uniqueness results based on ODE-leading-order terms 

Definition 2.2 of a quasilinear symmetric hyperbolic Fuchsian system, as well as the 
conditions which the singular initial value problem for such a system must satisfy if we 
wish to apply Theorem 2.4 and thereby guarantee the existence of a solution, involve the 
specified leading-order term uq just as crucially as they involve the exponent vector /i and 
the functions Si, £2 and N appearing in Eq. (2.1). In some applications, it is not easy 
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to determine which choices of uq (if any) lead to these conditions being satisfied. Here 
we discuss an approach which starts with the choice of a leading-order term of a very 
restricted type (which we label "ODE" -leading-order terms), and provides an alternate 
set of criteria for the existence of solutions to the singular initial value problem (with 
an ODE-leading-order term). This approach, presuming the criteria are satisfied, also 
systematically produces a sequence (possibly finite) of improved leading-order terms, 
which effectively serve as progressively higher order approximations to the solution of 
the singular initial value problem. We detail this approach here. 

As we see in Section 3.3.3 below, the ODE-leading-order term approach is very 
useful in our analysis of the T 2 -symmetric spacetimes. In particular, this approach 
plays a crucial role in our use of Fuchsian methods to obtain an optimal collection of 
T 2 -symmetric solutions of Einstein's equations with AVTD behavior. 



2.4.1 The ODE-Fuchsian operator and ODE-leading-order terms 

We start by defining the differential operator Lode(^o)["]> which plays a central role in 
carrying out this approach. Presuming that we are working with a specified quasilinear 
symmetric hyperbolic Fuchsian system Eq. (2.1) with specified (as yet arbitrary) leading- 
order term uq and with specified parameters 5, (i, qo and q, we define the ODE-Fuchsian 
operator as follows: 

L O de(u )[v] := Dv + S^{u )No(u )v. (2.42) 

Here we note that since (by Definition 2.2) S\ t o is invertible, it follows that £ode(^o)[ - ] 
is well-defined. We also note that, since £ode(wo)[-] does not involve any spatial deriva- 
tives, it is essentially a parametrized set of ordinary differential operators (one for each 
point x G T 1 ) rather than a partial differential operator (hence the "ODE" label). 

Although not necessary yet at this stage, we assume q > 3 and qo = q + 2 (consistent 
with Theorem 2.4) in all of what follows. In particular this guarantees that all maps 
(including their first spatial derivatives) are continuous with respect to x (as follows from 
the Sobolev embedding theorem). For example, the ODE operator above is well-defined 
at every spatial point x under this condition. 

We wish to write the Fuchsian PDE system Eq. (2.1) in terms of the operator 
Lode(uq)[-]. Recalling the operational form L (no + w)[uo + w] = F(uq)[w] for Eq. (2.1), 
and noting that we can relate the operators L(uo + w)[-] and £ode(^o)['] as follows 

L(u + w)[v] =Si(u + w)Loj)e(uo)[v] + S 2 (uq + w)td x v 

( _i _, \ ( 2 - 43 ) 

+ Si(u + w) [S 1 (u + w)N(u + w) - S 10 (u )N (u )j v, 

we find that we can write out Eq. (2.1) in the form 

Lobe(u )[u + w] = F de(uo)[w], (2.44) 
if we define the term on the right hand side of Eq. (2.44) as 

^ODe(«o)H :=S^ 1 (u + w)F(uo)[w] - S^iuo + w)S 2 (u + w)td x {u + w) 

\ ( 2 -45) 
S 1 (uq + w)N(u + w) - S l0 (u )N (u )) (uq + w). 
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The expression -Fode(^o)[ - ] is well-defined so long as 5 is sufficiently small so that 5i(uo + 
w) is invertible for any given w in B$^ q ^ s for some s > 0; we presume this is the case in 
all of what follows. 

As noted above, a key aspect of this approach is the selection of a special class of 
leading-order terms. 

Definition 2.15. A leading- order term uq is an ODE-leading- order term if it satis- 
fies the condition 

L O DE(u )[u ](t,x) = 0. (2.46) 

Expanding out the expression for £ode(^o)['L we see that an ODE-leading-order 
term no must satisfy Duo(t,x) + N(x)uo(t,x) = 0, where N := Siq(uq)Nq(uq) is (by 

definition) independent of t. For those very special cases in which TV is independent 
of uq, Eq. (2.46) is a parametrized set of linear ODEs, which can be readily solved 
for uq. More generally, Eq. (2.46) is nonlinear and therefore not so easy to analyze. 
We are interested here only in those cases in which we can establish that solutions to 
Eq. (2.46) exist. For those cases, we proceed to seek solutions of the singular initial 
value problem for the system Eq. (2.1) with ODE-leading-order term uq\ we call this the 
ODE-singular initial value problem. Observe that the solution to an ODE-singular 
initial value problem, if obtained, behaves in a way which suggests that as t — > 0, the 
spatial derivative terms in Eq. (2.1) become negligible. Indeed, this is true for the AVTD 
solutions of the Einstein field equations which we treat in Section 3 below. 

It is useful for our work below to notice that under the assumption Eq. (2.46), a 
combination of Eqs. (2.45) and (2.43) yields 

-Fode(«o)M = Si 1 ^ + w)$(u )[w] - S^iuo + w)S 2 (u + w)td x w 

/ , n (2.47) 

- iS 1 (u + w)N(u + w)- S lfi (uo)No(uo)J w, 

where 3~(«o)[w] is defined in Eq. (2.7). 



2.4.2 (Order n)-leading-order terms 

We now use the ODE-leading-order term uq (presuming that it exists) and the ODE- 
Fuchsian operator Lode(^o)["] to generate a (possibly finite) sequence (u n ) of "qualita- 
tive" solutions (in a sense described shortly) to the corresponding ODE-singular initial 
value problem; these (u n ) play an important role in establishing a set of conditions which 
are sufficient to show that this ODE-singular initial value problem does admit a solution 
(see Theorem 2.21 below). 

We first consider the x-parametrized set of linear inhomogeneous ODEs 

L ODE (u )[v](t,x) = f (t,x), (2.48) 

where uq is a fixed ODE-leading-order term and /o is a specified inhomogeneity (whose 
regularity we discuss below). If we use W(t,x) to denote a fundamental matrix for the 
linear homogeneous equation £ode(^o)M(*> x ) = 0) an d if we let (u* t i(x), . . . ,u* >c i(x)) 
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represent free data for the initial value problem at to £ (0,(5), then the general solution 
to Eq. (2.48) may be formally written as follows: 

v(t,x) = W(t,x){u* A {x),...,u*, n (x)) T + W{t,x) [ S - 1 W- 1 {s,x)f (s,x)ds. (2.49) 

Jo 

We may then formally define the operator 

M OBE (u )[f ](t,x):=W(t,x) f s- 1 W- 1 {s,x)f (s,x)ds, (2.50) 

Jo 

which, if it exists, maps a given source function /o to the particular solution w = 
Hode('Uo)[/o] of Eq. (2.48) determined by (m=„ i i(x), . . . , u* iC i(x)) = 0. We notice that the 
definition of this operator is invariant if the fundamental matrix W is replaced by an 
equivalent fundamental matrix W i— >■ W ■ M for any invertible d x ci-matrix M £ H qo . 

To proceed, we need to identify conditions which are sufficient for the existence of 
Hode («())[■]• Noting that we may always choose the free data (it #j i(x), . . . , u*^(x)) in 
such a way that the first term in Eq. (2.49) equals (the already specified) uq, we wish 
to also show that these same conditions are sufficient to guarantee that the second 
term in Eq. (2.49), i.e., Hode(wo)[/o]> is higher order in time as t approaches and 
therefore serves as a remainder term (in the sense of Definition 2.1) for the singular 
initial value problem. We state the needed conditions in the following lemma, which is 
readily checked. 

Lemma 2.16 (Existence and properties of Hode^oM'D- Suppose that a quasilinear 
symmetric hyperbolic Fuchsian system Eq. (2.1) has been chosen (Definition 2.2) for, in 
particular, a fixed sufficiently small parameter 5, for differentiability indices q > 3 and 
Qo = Q + 2, and for an ODE-leading- order term uq (Eq. (2.46) ). Suppose that S^qNq is 
of Jordan normal form. Then M.ode{uo)['] is well-defined on the domain ^s,u,q f or every 
smooth exponent vector v, so long as each component of v is strictly larger than the real 
part of the negative of the corresponding diagonal element (eigenvalue) of SiqNq. The 
target space of Mode(uq)[-] is X$^ q for any smooth exponent vector ju < u , and one has 
the estimate 

mODE(u )[fo}hfi, q < C5"\\f Q \\s,u, q , 

where C > depends only on the eigenvalues of S^^Nq, on the dimension d of the first- 
order system, on the choices of q, v, and on the difference between v and Jl; the constant 
p > only depends on the difference between v and Jl, and on q. 

In particular, we may choose Jl arbitrarily close to V, and we then check that the 
remainder w = Hode(uo)[/o], (the second term in Eq. (2.49)) is of higher order in t near 
t = (as measured by Jl) as the order of /o (measured by z7) becomes large. As well, as 
the difference between Jl and v diminishes, one may have to choose the constant C to 
be larger, and the constant p to be smaller. 

Some comments about this lemma are in order. For each quasilinear symmetric 
hyperbolic Fuchsian system and each choice of leading-order term uq, there exists an 
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invertible matrix T G H qo such that TS^^NqT^ is in Jordan normal form. We assume 
in the following that such a transformation has been applied to the system, and hence 
that S^qNq is in Jordan normal form. The fact that the matrices Si and 5*2 in the 
principal part are in general not symmetric after such a transformation has been carried 
out is not important for the arguments that follow. Moreover, for simplicity we assume 
for each exponent vector here that those of its components which correspond rto the 
same Jordan block of S^qNo have the same value. 

As a consequence of this result, we may formally define the following sequence: 

Definition 2.17 ((Order n)-leading-order sequence). Suppose q > 3. With wq = 0, we 
formally set 

w n := M-ODE(uo)[FoDE(uo)[w n -i]], (2-51) 
for all positive integers 1 < n < q — 2. The (order n) -leading- order terms are then 
defined by 

u n := u + w n , (2.52) 

for < n < q - 2. 

To turn this formal specification of the (order-n)-leading-order-sequence into a defi- 
nition, we need to state sufficient conditions for the composition on the right hand side 
of Eq. (2.51) to be well-defined for each n. We do this in the proposition below. This 
proposition also proves that the sequence is characterized by certain properties which 
are relevant to the two roles which it plays: i) an increasingly accurate sequence of 
approximations to the solution of the singular initial value problem with ODE-leading- 
order-term uq (presuming that such a solution exists) ; and ii) a sequence of "new" , and 
"better" leading-order terms which can be used to define new singular initial value prob- 
lems (closely tied to the original) for which we can prove the existence of solutions. The 
following proposition states the manner in which the first use makes sense, and provides 
the first step towards proving that the second use works. 

Proposition 2.18 (Existence and properties of the (order n)-leading-order terms u n ). 
Let q > 3 and qo = q + 2. Suppose that a quasilinear symmetric hyperbolic Fuchsian 
system Eq. (2.1) has been chosen satisfying Definition 2.2 for an ODE-leading- order 
term uq satisfying Eq. (2.46), for fixed parameters 5 (sufficiently small) and [i, and for 
all differentiability indices q' in the interval [3, q] . Here we require that the exponent 
matrix whose existence (in specifying the function spaces containing S± 1,^2,1, and 
N\) is a necessary part of the definition of a quasilinear symmetric hyperbolic Fuchsian 
system (see Definition 2.2), can be written as Cij = for some vector-valued exponent £ 
with strictly positive entries. Suppose that the matrix S^q(uq)Nq(uq) is given in Jordan 
normal form, and suppose in addition that the following conditions are satisfied for all 
5' G (0,5] and all integers q' G (3, q]: 

(i) The remainder exponent vector /i is strictly larger than the negative of the cor- 
responding diagonal elements (eigenvalues) of S((qNq and satisfies the modified 
block diagonality conditions: For every w G -B,5',/i,g',s; we have 

H\fi]Si(uQ + w) = Si(u + w)TZ[n\, K[n\N(u + w) = N(u + w)K[fi], 
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and, there exists r > 0, so that 

TZ[fi]tS 2 (u + w)K[-u] G B s , Xtq , >r , (2.53) 
/or £ defined in terms of £, as above. 

(ii) There exists an exponent vector u with u > /i and a constant r > 0, so i/iai ^(tto) 
maps B S i^ q i jS into B 5 i^ q ^ r . 

(Hi) For all w G Bg/^ q > s /2 and to G ^<s',/t,g',s/2 / or an 2/ exponent vector fl which satisfies 
p. > H and with respect to which the system is block diagonal, there exists a constant 
r > /or which 

2(uo)[w + w] - 3"(n )H G B S i^ +v _ m > ir , 
Si{u + w + uj)-S 1 (u + w)e B S /^_ M i )r , 
K[u]t(S 2 (u + io + w) - S 2 («o + G B s ,fi + £_^ q ^ r , 

N(u + w + u) - N(u + w) G Bv^+^^q,^. 

Moreover, there exists a constant C > suc/i £/ia£ t/ie norm of each of these 
quantities can be bounded as follows, 

\\3{U )[W + Lj] - 3(u ){w]\\ S >fi +u - M > < C\\u\\ S tfa q t, 

with analogous inequalities holding for Si, S2, and N. 

Then, the sequence u Tl specified in Definition 2.17 is well-defined and, for some 5 G (0, 5] 
and constants 7 > 0, one has 

u n -u e Bg^ q _ {n _ 1)}S/2 , u n - u n _i G % iH _ (n _i )7) ,_ (n _i W2 , (2-54) 

/or all 1 < n < q — 2. Moreover, the residual 1 of u n , defined by 

Res[u n ] := L(u n )[u n ] - f(u n ), (2.55) 

is contained in X$^ +nK ^ )q _ n . 

We make a few remarks here concerning some of the details of this proposition. First, 
we observe that as a consequence of the definition of this sequence, presuming that we 
start with an ODE-leading-order term uq of a certain order of differentiability, we find 
that the first term of the sequence U\ retains that regularity (up to order q), while the 
rest of the elements of the sequence (112,113...) generally do not. This is true because, 
since w = is smooth, it follows from the formula Eq. (2.47) that Fode(^o)[^o] has 
q derivatives; the same is then true for w\. But then since -FodeC^o) maps Xs^ )Q to 

1 The residual of a function u is defined so that u satisfies the system Eq. (2.1) if and only if its residual 
vanishes; i.e., Res[u] =0. So the residual measures the accuracy of the sequence (u n ) as approximate 
solutions to the system Eq. (2.1) 
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Xgp q _i, we see that u>2 has only q — 1 derivatives. The same loss of a derivative occurs 
for each successive element of the sequence. 

Secondly, we remark that the modified block diagonal conditions for /x are a slight 
generalization of Definition 2.3. In particular, it is not necessary that 1Z\p] commute 
with S2 here. We have chosen a formulation of Condition (i) which applies directly our 
applications — see Sections 2.4.4 and 3. Condition (i) can, however, be generalized to 
match more general situations. 

Thirdly, we note that Condition (iii) can be checked in our applications using the 
tools in the appendix; see Section B. 

Finally, we note that for the special choice of dj = £j used here, the space of matrix- 
valued functions Xg^ t q can be equivalently written as Xg t ^ q . This latter space — a Banach 
space of matrix-valued functions Xg^^ q with a vector-valued exponent £ — is defined in 
essentially the same way as above in Eq. (2.3); the key difference is that the norm used 
to define this space is the iif 9 -norm of - the matrix product of the matrix 7£[£] 

formed from £ (see Eq. (2.2)) times S. The motivation for specializing the matrix- valued 
exponent £ in this way and thence introducing the new notation Xg£ q for matrix-valued 
functions is to be able to express Condition (iii) of Proposition 2.18 in a natural way. 

The proof of Proposition 2.18 depends upon tight control of -FodeC^o)! - ] as given by 
Eq. (2.47), and tight control of the inverse of Sq. We obtain this needed control using 
the following two lemmas. 

Lemma 2.19. If the hypothesis for Proposition 2.18 holds (presuming as usual that 5 is 
sufficiently small so that S\(uq + w) is invertible for all w G Bg^qg and presuming that 
5' and q' satisfy the conditions stated in that hypothesis), it follows that for some r > 0, 
Si 1 (uq + w) G jB<s',o,g',r f or a M w £ Bfr ^ qi s . Moreover, the operator given by 

w 1 y 5f 1 (u + w) - (5i )0 (u )) _1 

maps Bg> u q' s into Bg'^^/.r for some constant r > 0, and this operator satisfies the 
difference condition 

Si 1 ^ + W+UJ) - S^Wo + W) G Bs'fi+^^q'^r 

for some exponent vector £ > 0, and for some constant C > it satisfies the inequality 

HSf^Uo + W + LO) - SfVo + W^U'fi+C-M' - C \\ U3 \\s',f£,q', 

for all w G Bgr a f 2 and all oj G Bgt^ q i >s i 2 ; here ju is any exponent vector which satisfies 
the inequality ju > {x, and for which the system is block diagonal. 

The proof of this lemma relies on i) the fact that the inversion of invertible matri- 
ces is a smooth map, ii) the fact that both TZ[p] and 1Z\fi] commute with S\, and iii) 
Proposition 3.9 in Chapter 13 of [40]; we omit the details here. 

Lemma 2.20. If the hypothesis for Proposition 2.18 holds, then there exist positive 
constants 7 and r so that for all constants 5' G (0, S] and all integers q' G [3, q\, and for 
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all exponent vectors ju > fx with respect to which the system is block diagonal, we have 
that F de(uq)[w] € B S r )fl+Jjq i^ ljr and 

F de{uq)[w + ui\- F O de(u )[w] G Byji+itf-if 

Further, there exists a constant C > such that for all w 6 ^s',ii,q',s/2 an d a H w G 
Bs',}i,q>,s/2, we have 

\\F O De{u )[w + Uj] - F De{uq)[w\\\s> fi^qi-l < C\\u)\\ S i firf . 

Proof of Lemma 2. 20. The first statement is easily obtained by multiplying the expres- 
sion Eq. (2.47) for -Fode('Wo)['] by the quantity S\(uq + w) and then using the facts that 
TZ[fj] commutes with Si(uq + w), and that Si(uo + w) is in -B^o^r for some r > 0, and 
also applying Lemma B.2 and Lemma 2.19. To prove the rest, we multiply this same 
expression for -Fode(^o)["] by 5i(«o + w), and then calculate 

Si(u + w)(F ODE (u )[w + uj] - Fode(^o)H) 
=3 f (u )[w + u]-3 f (u )[w] 

- (Si(u + w + uj) - S\{uq + w)) J(u )[w + uj] 

- (S2(uq + w + uj) - S 2 (u + w)) td x (w + uj) 

- S 2 {u + w)td x uj (2.56) 

- (N(uq + w + uj)- N(u + w)) (w + uj) 

- (N(u + w)-N (u ))uj 

- Si(u + w)(S^ 1 (u + w) - S^q(uo))N (u )uj 

- (S\(uo + w + uj) - Si(u + w)) N (u )(w + uj). 

Applying arguments of the sort used to verify Proposition 2.18 and Lemma 2.19 together 
with the estimates included in the hypothesis, we obtain the conclusion. □ 

We now proceed to prove Proposition 2.18: 

Proof of Proposition 2.18. We first show that the sequence iw n ) (and the corresponding 
sequence (u n )) is well-defined at least for finitely many sequence elements. It follows 
from Condition (ii) of the hypothesis that 3~(ito)[0] € B§ vqr for some r > 0. Noting 
(see Eq. (2.47)) that -Fode(^o) evaluated at reduces to 5^ 3"(uo)[0], we infer from 
Lemma 2.16 and Lemma B.2 that the term w\ is hence well-defined and is contained 
in Xg^q. It then follows from Lemma 2.20 and Lemma 2.16 (whose hypotheses are 
satisfied) that the operator Hode('Uo)[ - ] is well-defined, and consequently that w n is 
well-defined for all 2 < n < q — 2. These functions are all elements of -X^q-fn-i)- Using 
the estimate for the operator HIode(^o)['] stated in Lemma 2.16, we verify that if we 
shrink the time interval (0, 5] to (0, 5] as stated in the hypothesis of the Proposition under 
consideration, then we can show that finitely many of the sequence elements stay in a 
ball of fixed radius. We have thus verified the first statement appearing in Eq. (2.54). 
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Note that for convenience, in the remainder of this proof we continue to write 5 instead 
of 5; however, we reserve the right to repeatedly shrink the time interval as necessary (a 
finite number of times). 

We next argue by induction that the second statement in Eq. (2.54) holds. We 
presume that the differentiability index q is sufficiently large so that there exist nontrivial 
n < q— 2. To initialize the induction, we note that for n = 1, this statement says that u\ — 
u o G B &4iq s / 2 - Noting that, by definition, m - u = w\ = H de(^o)[-Fode(uo)[0]], and 
recalling from above that this term is contained in B^ ^ q s j2, we verify the initialization. 

To continue the induction argument, we suppose now that for some positive integer 
m < n there is an exponent vector /i(' m ) > fi such that w m — w m -i G B & q _r m _u s / 2 , 
and further suppose that the same is true for all positive integers i less than m (with cor- 
responding exponent vectors /u®). Using the definition of Wi, together with Lemma 2.16 

and Lemma 2.20, we find that there exists a pair of exponent vectors fi < ^^ m+1 ' < 
u (m+i) such that 

\\w m+ i - WmWs^im+V ,q-m = IPoDE^o) [-^ODE^o) [w m ] ~ F O Be(u ) [w m -i]] ||<5 )At (m+i) i9 _ m 

< C5 p \\Fo- DE (u )[w m \ - Fon^uo^Wm-ilWg^m+i)^^ 

< C5 p C\\w m - w m _i|| (5jjy ( m+ i)_ 7j(? _ (m _ 1) . 

In carrying out this calculation (with 7 being the quantity hypothesized in Lemma 2.20), 
we note that the operator HIode^o)! - ] is well-defined here according to Lemma 2.16 and 
Lemma 2.20 since w m — w ra _i G B s ^( m ) g _( m _i) s /2 with > fi. Finally we note that 
the constants C and p may depend in particular on q and m, but this dependence is not 
a problem for carrying out our argument since we are only interested in finitely many 
sequence elements. 

To complete the induction argument, we verify that since we have assumed (as part 
of the induction) that w m — w rn -i G ^^("^-(m-l)^^) it follows that so long as i/' m+1 ) — 
7 < /i*™' holds, we have the final right hand side of the above inequality finite. Therefore 
the initial left hand side must be finite, and this holds for any ^ m+l \ so long as //( m+1 ) < 
^(m) _j_ ry_ We satisfy these conditions by choosing ^( m+1 ) = p + rriK'y for any k < 1. 
Noting that this is the case for all m, with k chosen independently of m, we conclude 
that Eq. (2.54) holds, after having identified K7 with 7 to simplify the notation. 

It remains to verify that Eq. (2.55) holds for the residuals of the sequence (u n ). Using 
Eqs. (2.44), (2.46) and (2.51) we calculate 

Res[u n ] = L(u + w n )[u + w n ] - F(u )[w n ] 

= -Si(u + W n ) (-FoDE(^o)[^n] - F ODE (uo)[w n -i]) . 

Since w n — w n -\ G -B<5 iM +( n -i)7, (? -n+i,s/2! it follows from Lemma 2.20 that 

R,es[u n ] G Xg^j_ n ^^q_ n . 

□ 
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2.4.3 (Order n)-singular initial value problem 

Proposition 2.18 shows that, so long as we can find an ODE-leading-order term uo and 
so long as certain conditions hold, the difference u n +i — u n behaves like a power of t 
near t = 0, with this power increasing monotonically with n. It is hence meaningful to 
consider, in addition to the ODE-singular initial value problem with leading-order term 
uo, a sequence of (order n)-singular initial value problems which use (order n)-leading- 
order terms u n (n < q — 2). In view of the relationship between no and the sequence 
(u n ), we may write the same solution u of a given singular initial value problem either 
in the form u = uq + w for a remainder w in Xguq, or, as u = u n + u for a remainder 
uj in Xgti q with ju increasing suitably with n. The same can be done for any of the u m 
(m < n) in the (order n) -leading-order term sequence. 

We now use the (order n)-leading-order terms to argue that, at least for the smooth 
case (q = oo), if the conditions of Proposition 2.18 are met, the ODE-singular initial 
value problem, and correspondingly the (order n)-singular initial value problem have 
solutions. 

Theorem 2.21 (Existence and uniqueness for the ODE-singular initial value problem). 
Suppose that a quasilinear symmetric hyperbolic Fuchsian system with ODE-leading- 
order term uq has been chosen which satisfies the hypotheses of Proposition 2.18 for 
all finite values of (differentiability index) q. Then, for some sufficiently small 5\ £ 
(0,(5] and for a sufficiently large n, there exists a unique solution u of Eq. (2.1) with 
u — u n , D(u — u n ) 6 Xs 1:fJi + n ~f : oo> where u n is the (order n) -leading- order term defined in 
Definition 2.17 for this system. This solution u is also the only solution of the ODE- 
singular initial value problem with u — uq £ -Xfo n «>• 

This result states conditions which are sufficient for the ODE-singular initial value 
problem (with leading-order term uq) to admit (unique) solutions. In doing so, Theo- 
rem 2.21 provides a potentially very useful alternative to Theorem 2.4 of Section 2.2. 
Observe in particular that the hypothesis for Theorem 2.21 does not require that the 
energy dissipation matrix be positive definite with respect to fj,. 

Here we state and prove Theorem 2.21 only for the infinite differentiability case 
(q = oo). This smoothness restriction plays a role in the proof, since it allows one 
to always choose n large enough so that Condition (i) of Theorem 2.4 for the singular 
initial value problem with leading-order term u n is satisfied. If one tries to prove a result 
like Theorem 2.21 for finite differentiability order, then there is an upper bound for the 
possible choice of n, and consequently one may not be able to choose it large enough 
to satisfy Condition (i). However, in certain circumstances, a large but finite order of 
differentiability is in fact sufficient to carry through the proof. 

Proof. The basic idea of the proof is to reformulate the system using u n for the leading- 
order term in place of uo, and then verify that the hypothesis of Theorem 2.4 (in the 
case q = oo) is satisfied if n is chosen sufficiently large. To carry this through, we first 
argue that the system Eq. (2.1), which for the ODE-singular initial value problem can 
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be written as 

= L{u + w)[u +w]- F(u )[w], (2.57) 

can also be written as 

= L(u n + - ?(u n )[oj], (2.58) 

where we recall the definition Eq. (2.5) of the principal part operator L and the definition 
Eq. (2.7) for the operator 3~(u n )[-]. Here we use w for the remainder term corresponding 
to no and we use oj for the remainder term corresponding to u n (hence uq + w = u n + 
oj). To show this equivalence, we note the relations L(u n + oj)[v] = L(uq + w)[v] and 
F(u n )[uj] = F(uo)[w], and then using these we calculate 

= L(u + w) [uq + w]- F(u )[w] = L(u n + oj) [u n +u>]- F(u n ) [oj] 
= L(u n + oj) [u n ] + L(u n + oj) [oj] - F(u n ) [oj] 
= L(u n + oj)[oj] - 1{u n )[oj]. 

The equivalence of Eq. (2.57) and Eq. (2.58) immediately follows. 
We now choose a sequence of exponent vectors //W which satisfy 

H < < n + (n - 1) 7 , (2.59) 

and which are consistent with the block diagonal condition for Eq. (2.58); we note that 
this is possible for all sufficiently large integers n. Examining the singular initial value 
problem corresponding to Eq. (2.58), we verify that for any given sufficiently large integer 
n, this PDE system, together with u n as leading order term and exponent vector J^ n \ 
satisfies the conditions to be a quasilinear symmetric hyperbolic Fuchsian system. We 
also verify, based on Eq. (2.6), that for sufficiently large n (and therefore sufficiently 
large [i + (n — 1)7) the exponent vectors Jl^ can be chosen large enough to guarantee 
that the energy dissipation matrix Mq is positive definite. Consequently, this system 
satisfies Condition (i) of the hypothesis of Theorem 2.4. 

To check that Conditions (ii) and (iii) of Theorem 2.4 are also satisfied, we examine 
the operator 3~(u n ). Using Eqs. (2.43) and (2.47) together with Definition 2.17, we 
calculate 

J(u n )[oj] =S 1 (u + w) (F ODE (u )[w] - FoDE(«o)[w>n-i]) + S 2 (u + w)td x oj 
+ 5i(«o + w) (sr Vo + w)(N(u + w)- N (u )) 

+ {S^(u + w)- S^(uo))N (uo)^oj. 

By comparing the first line of this expression with Eq. (2.56), we notice that all 
spatial derivative terms cancel; hence there is no loss of regularity in this expression as 
is the case for the operator -Fode(^o) itself. We therefore get estimates analogous to 
those in Lemma 2.20, with q — 1 replaced by q. 

Combining the assumptions for Si, S2 and TV" which are stated in the hypothesis 
of Proposition 2.18 (and therefore included in the hypothesis of Theorem 2.21) with 
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the upper bound stated in Eq. (2.59), we readily show that all of the conditions of 
Theorem 2.4 hold for q = oo. Then the consequent application of Theorem 2.4 shows 
that so long as n is sufficiently large, there exists exactly one solution u = u n + uj with 

U G X S! ,0O" 

We wish to show next that for such a fixed chosen value of n, in fact uj G Xg lt u+nj,oo- 
To show this, we consider an integer n + which is large enough so that Jl^ n+ ^ > /i+(n— 1)7. 
Applying the same argument as above, but now with n + instead of n (and hence using 
u n+ as the leading-order term), we obtain a solution u = u n+ + uj which has the property 
that uj G ~(n + ) OQ . Uniqueness of the singular initial value problem with respect to 
u n implies that u equals u. Moreover, we have 

w = w n+ -w n + oj = (w n+ i - w n ) + ... + (w n+ - W n+ -l) + uo. 

Given that w n+ i — w n S Xg lt ^ +nr y !00 , we obtain the desired result 

To conclude the proof of this theorem, we must show that any solution u of the form 
u = uq + w with w G Xg 1 n,oo must equal the solution u. To show this, it is useful to write 
u = u n + w — w n , where u n is defined by Eq. (2.52) and w n is defined by Eq. (2.51). Then 
if we can verify that w — w n G ^SxM n ),ooJ ^ follows from uniqueness that uj = w — w n , 
and hence that u = u. We make this verification by using induction to show that, in 
fact, w — w m G Xj li/i+m7]00 holds for every non-negative integer m. In the case m = 0, 
we have w — u>o = w G Xg 1 ^ iOC which implies the claim for m = 0. Suppose the claim 
has been shown for m = niQ > 1. We know that w is a solution of the equation 

-^ode^o)^] = Fode(u )[w], 

while w mo+ i is a solution of 

LoDE(uo)[w mo+ i] = ^ODE(«o)[^m ]- 
Taking the difference, we obtain 

L ODE (u )[w - w mo+ i] = F ODE (uo)[w} - F DE(uo)[w mo }. 
We can write this formally as 

w - w mo+ i = Mode(uo)[F OD e(uo)[w] - FoDE(^o)[^mo]]- 

Now, the fact that w-w niQ G Xs ufl+m07:OO implies that F -de(uo)[w]- F ODE(uo)[wm ] G 

At +( m0 +i) 7jOO (Lemma 2.20). Consequently (see Lemma 2.16), the operator Hode(wo)["] 
is well-defined. This completes the proof. □ 
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2.4.4 An example: the Euler Poisson Darboux equation 

We consider now the example of the Euler-Poisson-Darboux equation (see also [2] for 
another example) 

D 2 u(t, x) - t 2 u xx (t, x) = f (t, x). (2.60) 

Here, u(t, x) is the unknown (assumed to be a scalar function), and fo(t, x) is a specified 
scalar function. The Euler-Poisson-Darboux equation is second order, and in previous 
work by two of the authors [13] on semilinear second-order Fuchsian systems, it has been 
shown that this equation admits unique solutions to the singular initial value problem 
with leading-order term 

u (t, x) = u*(x) log t + u**(x), (2-61) 

(for arbitrary functions i** and u**) so long as /o = 0(t v ) with v > 0. We seek to show 
that we obtain these same results using the first-order methods which we have developed 
here. In particular, this example demonstrates the usefulness of the techniques discussed 
in Section 2.4, thereby serving as a linear warmup example with which we can explore 
some of the issues which arise below in our discussion of the application of these methods 
to the fully nonlinear T 2 -symmetric Einstein's vacuum equations in Section 3. 

To apply the first-order theory developed in this paper, we first convert this equation 
into a first-order system by setting 

u± := u, U2 := Du, u%:=td x u, U := (i*i, 1*2, Us) T . (2.62) 

Eq. (2.60) then takes the form of a first-order evolution system 

SiDU + S 2 td x U + NU = /, (2.63) 

with 

/0 \ 
5i = Diag(l,l,l), S 2 = -1 , N = 

\0 -1 oj 

plus a constraint equation 

A u := u 3 /t - 8 xUl = 0. (2.64) 

Observe that in working with the Euler-Poisson-Darboux system in this first-order form, 
one first treats the components 1*1,1*2, and 1*3 as independent functions whose evolution 
is determined by Eq. (2.63). This means that we solve the singular initial value problem 
of this system with respect to a leading-order term motivated by Eq. (2.61). Then, in a 
second step, we identify 1*1 with the original variable it and consider the two remaining 
relations Eq. (2.62) as constraints: the one involving the time derivative is automatically 
implied by the first of Eqs. (2.63) (the evolution equation for 1*1), while the one involving 
the spatial derivative gives rise to the condition A u = in Eq. (2.64). Let us start with 
the first step. 

One readily verifies that this evolution system is of (quasilinear) symmetric hyper- 
bolic Fuchsian form for any choice of leading-order term, and hence our theory can, in 
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principle, be applied. Our approach is to find a leading-order term for the first-order 
variables which is consistent with Eq. (2.61) and which, in addition, is an ODE-leading- 
order term. We easily determine that the general solution to Eq. (2.46) for Eq. (2.63) 
takes the form 

Uq = (Ci + C 2 logt, C 2 , C 3 tf, (2.65) 

for the spatially-dependent parameters C\(x), C 2 (x) and C%{x). However, we see that 
this leading-order term can only be consistent with Eqs. (2.61) and (2.62) in the special 
case k* = and C 2 = 0. Hence, this approach for finding a leading-order term fails. 

We circumvent this problem as follows. For a specified function u* (which is at least 
second order differentiable; we specify its necessary regularity more precisely below), we 
define 

u := u — u*{x) logt, (2.66) 

and work with the evolution equation for u rather than that for u. Substituting Eq. (2.66) 
into Eq. (2.60), we obtain 

D 2 u - t 2 u xx = t 2 login" + f (t, x), 

where u'l indicates the second derivative of . Now, setting 

Si := u, u 2 := Du, u 3 := td x u, U := (ui,u 2 ,u 3 ) T , 

we obtain the evolution equation 

S^U + S 2 td x U + NU = /, (2.67) 

for the same matrices Si, <S 2 , ./V as above, but with 

/ = (0,/ + t 2 logt<,0) T 

In terms of u, the constraint Eq. (2.64) takes the form 

A s := %/t - d x ui = 0. (2.68) 

Choosing the ODE-leading-order term for the U formulation to be of the same form as 
Eq. (2.65), we have 

U = (C 1 + C 2 logt,C 2 ,C 3 t) T , 

but now (in view of Eq. (2.61)) we are led to choose the parameter functions in the form 
C\ = u** , C 2 = and C 3 = u' m ; hence 

U (t,x) = (u**(x),(W„(:r)) T (2.69) 

The function appearing in Eq. (2.66) together with the function u** introduced here 
together comprise the full range of free data suggested by Eq. (2.61). Both play the role 
of asymptotic data functions. 
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Having found a suitable representation of the equations and the leading-order term, 
we write the unknown U of the evolution system as U = Uq + W, and look for sufficient 
conditions for the existence of solutions to the singular initial value problem in this form, 
with W as a remainder term. To enforce the remainder falloff properties, we choose an 
exponent vector \i = (/^i, fj,2, M3) and, in view of Eq. (2.69), we require that 1x1,^2 > 
and /i3 > 1. 

We first seek to prove existence of solutions using Theorem 2.4. To satisfy the 
block diagonality condition of Theorem 2.4 we must set \i\ = 1x2 = /^3- We therefore 
simplify the notation by writing the exponent vector as (//, //, /i) for some smooth scalar 
function \i which, from above considerations, must be greater than one. Observe here 
that, while this equality of all components of the exponent vector is necessary to satisfy 
the hypothesis of Theorem 2.4, it does appear to be an artificial restriction. Under 
reasonable regularity assumptions, we might rather expect that if the first and second 
components are 0{t^), then the third component of W should be 0(t^ +1 logt); the logi 
factor may arise from derivatives of since fx is generally not constant. In any case, we 
readily verify that the energy dissipation matrix 

(H -1 \ 
M = ix , 
\0 fx-lj 

is positive definite so long as fx > 1. 
Calculating 

nU )[W] = F(U )[W] - L(U + W)[U ] = (0, /„ + i 2 (logt < + <J, 0) T , 

we now suppose that W G Xs,(^,^,^),q an d /o G Xg^ q for v > 1. Then 3~(J7o)[^] £ 
Xs,(v,i/,v),q ^ G -ff 9+2 (r 1 ), where j/ = i/, if i/ < 2, or, we have < 2, if v > 2. 

Choosing q > 3, we verify that Theorem 2.4 implies the existence of solutions of the 
evolution system U = Uq + W with W G Xg 1 j /J/ ^ f j / \ (j for <5i sufficiently small 1 and for an 
exponent n G (l,min{2, u}). For any specified set of the asymptotic data u* and u**, 
we find that the solution is unique for remainders in the space X^^^^^. 

Given any such solution of the first-order evolution system, our next step is to identify 
u\ with u—\og(t)u* and then, if the remaining constraint = is satisfied, to conclude 
that u is actually a solution of the original second-order equation Eq. (2.60) with leading- 
order term uq = u* log t + u** and with remainder w = w\ (the first component of the 
vector W) in X& 1% ^ q . To determine if the constraint is satisfied, we use the evolution 
equation Eq. (2.67) to calculate the time derivative of the constraint violation term A®, 
obtaining 

DA 2 = 0. (2.70) 

We then note that (i) if we construct A^ using Uq from Eq. (2.69), we get A^ q = 0; 
and (ii) if we combine the evolution equation Eq. (2.70) with the leading order term 

1 In fact, since the equations are linear, we can extend the solution to all positive times t > 0. 
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A^ as well as q > 3 and other appropriate choices of /i, etc., then we find that Ag 
satisfies a singular initial value problem which satisfies the hypothesis of Theorem 2.4. 
Noting that A^ = is a solution to this singular initial value problem, and recalling 
that Theorem 2.4 implies that solutions are unique, we see that indeed, the constraint 
Ag = must be satisfied. 

While this approach to analyzing the singular initial value problem for the Euler- 
Poisson-Darboux system does produce a solution, it is unsatisfactory for two reasons. 
First, it does not allow us to treat the case in which /q € ^s,u,q f° r v < \. Second, if 
v > 1, this approach does not exclude the possible existence of other solutions u with 
remainders w in Xg tfltq for jjl < 1. Both of these issues are resolved if we use an alternative 
approach based on Theorem 2.21 and the use of (order n)-leading order terms. In doing 
this, we pay a price in that we must require a that the spatial derivative parameter q is 
infinite. 

If we are to work with Theorem 2.21, a key requirement is that we start with an ODE- 
leading-order term; we have already fulfilled this requirement by our choice of Uq. We now 
have the advantage that we do not need to impose the block diagonal condition, but only 
the modified block diagonal conditions, see Condition (i) in Proposition 2.18, and also not 
the positivity of the energy dissipation matrix in choosing the remainder exponent vector 
fi; we may work with \i = (m, (J,2, ^3) for any fi±, ^2 > and 1 < /X3 < [12 + 1, thereby 
permitting the full range of values of fj, for which the singular initial value problem is 
meaningful. Notice that the upper bound for /X3 is implied by Eq. (2.53) and is related 
to the observation above that a spatial derivative of a spatially dependent power of t 
may introduce additional logt-terms. Proceeding, we suppose that we have chosen some 
fo E Xsfi^ with v > 0. Any choice of fi satisfying the above conditions is consistent 
with Condition (ii) of Proposition 2.18 (as part of Theorem 2.21) if ^2 < min{2,P'}. 
Choosing it*, u** G C 00 ^ 1 ), we then verify straightforwardly that Condition (iii) of 
Proposition 2.18 is satisfied. We conclude that there exists a solution U of the evolution 
system with U—U n £ ^<5 1 ,( Atl)M2 ,M3)+n7,oo ^ or some constants <5i > and a sufficiently large 
integer n. This solution is unique, with the remainder U—Uo contained in ^ ( 5 l! ( /il)At2ijU3 ) ! oo- 

Having verified the existence of solutions to the first-order evolution system, we wish 
to show again that the corresponding solution is actually a solution of the original second- 
order equation by considering the constraint Eq. (2.68). This can be done essentially as 
discussed above. 

To illustrate the use of the leading-order term approach, we choose the source term 
in the form fo(t,x) = f*{x)t 1 / 2 for a smooth function /*, and calculate 



/4/„(x)< 
U 1 = U Q +\ 2U{x)t l 




and 





U 2 = U l + 







\^4t 3 / 2 i gt/: + ^ 3 (io g ^i 3) + 



u 
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One may continue to calculate the sequence, and one verifies (in accord with the last 
statement in Proposition 2.18) that the residuals corresponding to this sequence are 
contained in X spaces of monotonically increasing exponent. 

3 T 2 — symmetric vacuum Einstein spacetimes 

3.1 Objective of this section 

As noted in the Introduction, one of the main motivations for this work is to explore the 
singular regions of certain classes of solutions of the Einstein gravitational field equa- 
tions. In particular, as a step towards studying the strong cosmic censorship conjecture 
in families of solutions characterized by relatively large isometry groups, we use the 
Fuchsian formulations developed here to show that there are large sets of solutions in 
these families which exhibit AVTD behavior in a neighborhood of their singularity. 

We work here with spacetimes which are characterized by a spatially-acting T 2 isom- 
etry group, but do not have the further restriction of a non- vanishing "twist", which 
defines the familiar Gowdy spacetimes. Following convention, we refer to them as the 
"T 2 -symmetric spacetimes"; if they also satisfy the Einstein equations, we call them 
"T 2 -symmetric solutions". While much is known regarding the Gowdy spacetimes, in- 
cluding a proof that strong cosmic censorship holds for the Gowdy spacetimes with T 3 
spatial topology [35] and for polarized Gowdy spacetimes with any allowed spatial topol- 
ogy [18], much less is known about the T 2 -symmetric solutions. For both the Gowdy 
and T 2 -symmetric families, the presence of the T 2 isometry effectively reduces the anal- 
ysis to that of a PDE system on a 1 + 1 dimensional manifold. One notable difference, 
however, is that while the Gowdy PDE system is semilinear, that of the T 2 -symmetric 
solutions is quasilinear. 

The first work showing that there are (non-polarized) Gowdy spacetimes with AVTD 
behavior is that of Kichenassamy and Rendall [30] which uses Fuchsian methods to show 
that this is true for analytic Gowdy solutions on T 3 . The later work of Rendall [34] shows 
this for Gowdy spacetimes which are smooth, again using Fuchsian methods (adapted to 
smooth solutions rather than analytic solutions). Fuchsian methods have been used [24] 
to verify that there are analytic polarized T 2 -symmetric solutions with AVTD behavior. 
Here, we use the results presented above to show the same for T 2 — symmetric solutions 
(polarized and half-polarized) which are not analytic. 

3.2 T 2 -symmetric spacetimes 

The family of vacuum T 2 -symmetric spacetimes is characterized by a T 2 isometry group 
which acts effectively on each spacetime in the family, with the generating Killing vector 
fields being everywhere spacelike. We assume that each such spacetime is the maximal 
globally hyperbolic development of an initial data set on a compact Cauchy surface, with 
the data invariant under an effective T 2 action. One more condition distinguishes the 
spacetimes we consider here from the Gowdy subfamily. Let Y and Z be the generators 
of the T 2 isometry. The Gowdy subfamily is characterized by the assumption that 
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the distribution denned by the tangent planes orthogonal to the generators Y and Z is 
integrable. This condition is usually expressed as the vanishing of the two twists Ky and 
Kz- We work here with T 2 -spacetimes with at least one non- vanishing twist. Chrusciel 
has shown [17] that the vacuum Einstein equations force the twists to be constants, 
and that the condition of non-vanishing twist implies the Cauchy surfaces must have T 3 
topology. 

Such spacetimes can be foliated by areal coordinates, in which the time coordinate 
labeling each symmetry group orbit is equal to the area of that orbit. This coordinate 
system conveniently locates the singularity at t = except in the special case of flat 
Kasner, as is shown by Isenberg and Weaver in [27]. Local existence for these coordinates 
is shown by Chrusciel, [17], and global existence is proved by Berger et. al. in [8], and 
further clarified in [27]. 

Let y, z be coordinates on T 2 , and let x be the remaining spatial coordinate, which 
takes values in S l . The metric can be written [17] in the form 2 

g = e 2 ^- u )[-adt 2 + dx 2 ^+e 2U (dy + Adz+(Gi + AG 2 ^dx^ + e~ 2L T t 2 '(dz + G 2 dxj , 

where all the metric functions {?/, U, a, A, G\, G 2 } depend only on t and x. 

If both twist constants vanish, then the function a can be chosen to be a constant, 
in which case the above metric reduces to the Gowdy metric [22]. 

The polarized class of T 2 -symmetric spacetimes results from setting A equal to a 
constant in the initial data (or, equivalently, assuming that the dot product of the gen- 
erators Y, Z is initially the same at all spatial points 3 ), and verifying that this condition 
is preserved under evolution. While the polarized spacetimes are characterized by a ge- 
ometric condition, another subclass we consider, called the half-polarized T 2 -symmetric 
spacetimes, is defined by a restriction on the asymptotic behavior of the fields (see 
Section 3.3.1). 

Before writing down the Einstein vacuum equations, we make a few further coor- 
dinate choices to simplify the presentation. Without loss of generality we choose the 
generators such that Ky = 0, Kz = K 0. This can be achieved by choosing an appro- 
priate linear combination of any generators for the T 2 action. It is sufficient to consider 
K > since the transformation K — > —K preserves all conditions imposed thus far. 
Next we choose coordinates y, z on T 2 so that Y = d y and Z = d z . This can be done 
without changing the form of the metric above. Implementing these simplifications, and 
using the short-hand notation Ut '■= dtU for derivatives, we write the Einstein equations 

Let £ : = g(Y, •), C : = g{Z, •) be the generating forms of the distribution D. Frobenius' theorem states 
that D is integrable if and only if Ky '■= *d£ A £ A ( and Kz := *dC A £; A ( both vanish. 

2 In some representations of these metrics, two "shift constants" M x and M v also appear. These can 
be removed using gauge choices; hence we leave them out. 

3 Observe that if the dot product is constant, one can always find a new pair of generators (by taking 
linear combinations) which are orthogonal. 
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as the following system of PDEs, which includes a set of second order equations 



U tt + — -aU xx = -^- + ^- + ^2V A t - aA x), (3-1) 

= ^ + ^-4A t U t + 4aA x U x , (3.2) 
i 2 2a 

a x Vx . am a 2 x a xx 2 2 

»7tt - arj xx = — ^— + — - ^ + -y- - I7 t + al/ XJ (3.3) 



a set of first-order equations 



P 4C/ „2t/ 



7? t = tUf + taU 2 x +—{A 2 t +aA 2 x ) + ^aK\ (3.4) 
e 4C/ a 

rj x = 2tU t U x + —A t A x --^, (3.5) 
2t 2a 

e 2r > 

a t = --^a 2 K 2 , (3.6) 
plus a set of auxiliary equations 

G u = e 2ri ^a~AKt- 3 , G 2t = -e 2ri V^Kt~ 3 . (3.7) 

Here, the auxiliary equations originate from the definition of the twist constants Ky and 
Kz and from the "gauge" simplification Ky = noted above. 

Observe that the T 2 -symmetric Einstein system reduces to the Gowdy system in 
the standard areal coordinates if we set K = 0, a = 1, G\ = 0, and G% = 0. The 
Einstein equations in the Gowdy class are semilinear and a Fuchsian analysis with ana- 
lytic asymptotic data has been carried out by Kichenassamy and Rendall [30], and with 
smooth asymptotic data by Rendall [34] and by Beyer and LeFloch [12]. 

3.3 Existence of AVTD solutions to the Einstein vacuum equations 
3.3.1 AVTD behavior and heuristics 

What is the behavior of a singular solution to Einstein's equations near the singularity? 
In principle the behavior could be very complicated for a solution to a system of nonlinear 
PDE such as the Einstein equations. In [32, 6, 7] Belinskii, Khalatnikov, and Lifshitz 
(BKL) propose that generically the spacetime dynamics near the singularity is vacuum 
dominated, local, and oscillatory. According to this picture, an observer traveling toward 
the singularity (either backward or forward in time, depending upon the location of the 
singularity) would experience an infinite sequence of Kasner epochs, and each observer 
at different spatial points would experience a different, generally unrelated, sequence. 

Numerical simulations of T 2 -symmetric spacetimes [5, 9, 33] support this picture, 
except perhaps at points where spikes occur. Whether the complicated behavior found 
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near spikes, and the apparent prevalence of spikes, invalidates the BKL picture for 
general T 2 -symmetric solutions is far from clear. However, for the restricted family of 
polarized T 2 -symmetric solutions, numerical simulations indicate that a special form 
of BKL behavior occurs near singularities-asymptotically velocity term dominated, or 
AVTD, behavior-which is not dominated by what happens near spikes. In a spacetime 
with AVTD behavior, each observer experiences only a finite sequence of Kasner epochs 
in the approach to the singularity [25, 21, 23], and the limiting spacetime is different for 
each observer. 

While there are no analytical studies of inhomogeneous cosmological solutions which 
either confirm or deny the presence of general BKL behavior, as noted above there has 
been a significant amount of such work supporting the generic presence of AVTD be- 
havior in restricted families of solutions. Studies based on singular initial value problem 
formulations of Fuchsian PDEs are particularly well-adapted to doing this, since they 
involve specifying a choice of asymptotic behavior (a Kasner evolution independently at 
each point), and showing that there are solutions of the equations which approach this 
asymptotic behavior. If we can show that the Einstein equations for the polarized T 2 - 
symmetric spacetimes, together with certain choices of the leading order term, satisfy 
the conditions of the hypothesis of either Theorem 2.4 or Theorem 2.21, then we have 
confirmation that there are such spacetimes which have AVTD behavior. 

Observe that finding solutions in a given family of spacetimes with AVTD behavior 
does not imply that there are not solutions in that same family with a very different form 
of asymptotic behavior. However, since numerical simulations support AVTD behavior 
being generic among polarized T 2 -symmetric solutions, there have been no searches for 
alternative forms of asymptotic behavior among them. 

The name "asymptotically velocity term dominated" refers to the fact that the lead- 
ing order terms are chosen as asymptotic solutions of the "velocity term dominated" 
(VTD) system, which is formed from the Einstein equations by dropping terms with 
spatial derivatives. This step encodes the local aspect of the BKL proposal. It can be 
shown [24, 20] that the following expansions for the metric functions below asymptoti- 
cally solve this VTD system in the limit t — >■ 0. We write these expansions in terms of 
asymptotic data {k, £/**, A*, A** , 77* , a* , Gu, G2*} with the regularity of the data speci- 
fied below. 



U(t,x) 


= 2 (1 ~ k(x))log 


t + U**(x) + . . . , 


(3.8) 


A(t,x) 


= ^.^(x) -\- A** (x) 


t 2k ^ + ..., 


(3.9) 


7](t,x) 


= ±(l-k(x)) 2 \o. 


gt + rj*(x) + . . . , 


(3.10) 


a(t, x) 


= a*{x) + . . . , 




(3.11) 


Gi(t,x) 


= Gu(x) + ... , 




(3.12) 


G 2 {t,x) 


= G 2 *(x) + .... 




(3.13) 



Of particular importance here is the function k. It determines the Kasner exponents 
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Pi,P2,P3 °f the local Kasner solutions which are approached at any spatial point 
Pi = (k 2 -l)/(k 2 + 3), p 2 = 2(1 - k)/(k 2 + 3), p3 = 2(l + *)/(fc 2 + 3). 

We recall here that a T 2 -symmetric solution is defined to be polarized if the two 
Killing vectors corresponding to the T 2 isometry can be chosen to be orthogonal every- 
where. This is the case if and only if the metric coefficient A = const. A solution with 
AVTD behavior has this property if and only if the asymptotic data corresponding to 
A satisfy the conditions A** = and A* = const. Since A* = const can be gauged to 
A* = 0, we see that in the polarized case, there is effectively no free asymptotic data to 
choose which relates to A. There is an interesting relationship between the polarization 
condition and the sign of k: Examining equations (3.8)-(3.13), we find that if a solution 
is not polarized and has AVTD behavior, then there is power law blow-up at the singu- 
larity if and only if k is negative. Yet if that spacetime is polarized, then regardless of 
the sign of k, there is no power law blow-up at the singularity. 

The polarization condition is relevant to our application of our Fuchsian results to 
T 2 — symmetric solutions since, as we see below, our results cannot be applied unless the 
condition d x A* = holds for the asymptotic data. For polarized T 2 -symmetric solutions, 
this restriction on A* is automatic. It is important to note, however, that requiring 
d x A* = does not restrict us to polarized solutions. We may consider asymptotic data 
which has this restriction on A* , but has no restriction on ^4** . T 2 -symmetric solutions 
which are AVTD and which have asymptotic data of this sort are known to exist, and 
have been called "half-polarized" 1 [20]. Extending the results of both [24] (analytic and 
polarized) and [20] (higher regularity), we show here that there are large families of both 
half-polarized and polarized T 2 -symmetric solutions which are smooth or of even lower 
regularity, and which have AVTD behavior near their cosmological singularities. 

A general (neither polarized nor half-polarized) T 2 -symmetric solution, were it to 
be AVTD, would have asymptotic data with both A* and A** non-vanishing and non- 
constant. Based on numerical and heuristic considerations, however, it is expected that 
spacetimes with non-constant A* do not generally show AVTD behavior. Rather, these 
are expected to show Mixmaster-like BKL behavior at the t = singularity, or behavior 
which is even more complicated (with strong spike influence). We do not address this 
issue here. 

We now discuss two applications of our Fuchsian results which verify AVTD behavior 
in T 2 -symmetric solutions. For the first one, Theorem 3.1, we make only minimal 
assumptions regarding the regularity of the asymptotic data. The price to pay for this 
is that the result does not cover the full expected range for the function k = k{x) in 
Eqs. (3.8) - (3.13). For the second result Section 3.3.3, Theorem 3.6, we add regularity 
restrictions, but we do get the expected full range of allowed values for k. 

lr The use of the term "half-polarized" to describe AVTD solutions with "half of the asymptotic data 
of one of the gravitational degrees of freedom turned off first appears in a discussion of [/(l)-symmetric 
solutions with AVTD behavior, in [26]. 
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3.3.2 Existence of low regularity solutions with AVTD behavior 

The low regularity result, which we formulate, discuss, and prove in this subsection, 
is an application of Theorem 2.4 to the polarized and half-polarized solutions of the 
T 2 -symmetric equations. 

Theorem 3.1 (First result: AVTD (half)-polarized T 2 -symmetric vacuum solutions - 
finite differentiability). Suppose one chooses a twist constant K 6t, a pair of asymptotic 
data constants A* and rjo, and a set of asymptotic data functions k, ?7**,a* G i? 9+2 (T 1 ) 
(with a*(x) > 0), A** G H q+l (T l ) and Gu,G 2 * G H q {T 1 ) for any q > 3, which satisfy 
the integrability condition 1 



'"((1- k(x))UUx) - l(loga,)'(aO)ds = 0, 



together with, at each point x G T , either 

(i) k(x) > 1 + \/6 for arbitrary A** (the half-polarized case), 

(ii) k(x) > 1 + y/E or k{x) < 1 — v6 for A** = (the polarized case). 

Then there exists a 5 > 0, and a T 2 -symmetric solution U , A, rj, a, G\, G2 of Einstein's 
vacuum field equations with twist K of the form 

(U,A,Tj,a,G 1 ,G 2 ) = (U ,A ,r)o,a ,Gi :0 ,G2,o) + W. 
Here, the leading-order term (Uq, Aq, t]q, ao, Gi,o, ^2,0) is given by Eqs. (3.8)-(3.13), with 



r)*(x) := rjo + 



£ ((1 - k(X))UUX) - i(loga*)'(X)) dX. (3.14) 



The remainder W is contained in Xg (and DW G Xg a,q—i) for any exponent vector 
A* = (m, H2, V>z, Ha, Hh, V*) with 

1 < m(x) < min{2, (k(x) - 3)(k(x) + l)/2}, 



(2k(x) + yJl + Ak{xY)/2 < fi 2 (x) < 1 + 2k(x), 

< fj, 3 (x) < m(x), 

< m(x), fis(x),/i6(x) < (k(x) - 3)(k(x) + l)/2. 

This solution is unique among all solutions with the same leading-order term Uq and 
with remainder W G Xg^^ q . 



1 This integrability condition results from integrating the constraint equation Eq. (3.5) along a constant 
t circle, and taking the limit as t — > 0. Observe that here and below we use the notation U' tt ,(x), (log a*)', 
etc., to denote derivatives of functions which depend on x only. 
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Observe that by taking time derivatives of the Einstein field equations, we can also 
obtain corresponding statements about the behavior of a certain number of time deriva- 
tives D m W of the remainder function W. We do not elaborate on this any further 
here. 

This result, based on Theorem 2.4, does not imply uniqueness of the solutions within 
the whole class of solutions of interest: For a given choice of asymptotic data, Theo- 
rem 3.1 determines that there is exactly one solution with remainder W in spaces Xgn q 
with jj, given by Eq. (3.15). The full class of remainders compatible with the leading- 
order behavior Eqs. (3.8)-(3.13) however corresponds to exponents 

(j,!, (J, 2 - 2/c, /U 3 , /i 4 , ns,V>6 > o. 

Hence, for given asymptotic data there may exist further solutions in such a larger space. 
Strict uniqueness can be explored further using techniques involving (order n)-leading 
order terms. We return to this issue in Section 3.3.3 below; the price which we have 
to pay for strict uniqueness is that we need to require higher differentiability for the 
asymptotic data. 

In proving Theorem 3.1, it is useful to arrange the T 2 -symmetric Einstein vacuum 
equations, Eqs. (3.1)-(3.7), as well as the field variables, in a certain hierarchical form: 
Eqs. (3.1), (3.2), (3.4) and (3.6) together form a coupled evolution system (which we 
label the "main evolution equations") for the variables U,A,rj, and a. Eq. (3.5) serves 
as a constraint equation for this system, while Eq. (3.3) is effectively redundant, and 
can be ignored. The remaining equations Eqs. (3.7) are evolution equations for G\ and 
C?2, and can be handled after the analysis of the main evolution equations. 

We proceed now to focus on the main evolution equations, with the primary existence 
result for them - the main step toward a proof of Theorem 3.1 - being Proposition 3.2. 

Main evolution equations. To rewrite the main evolution equations as a first order 
symmetric hyperbolic Fuchsian system, it is useful to define certain new variables. Some 
of the choices of these variables are motivated by considerations in [24], others by the 
discussion above in Section 2.4.4. First, we set 

£ := d x a, (3.16) 

whose evolution equation is obtained by taking the spatial derivative of Eq. (3.6) and 
by substituting any occurrence of r] x by the constraint Eq. (3.5). One obtains 

e 2r? 

6 = —^ aR2 (*£ + a ( e A ^A t + 4t 2 U x U t )) . 

In all other evolution equations we use Eq. (3.6) to eliminate at and replace a x by £. 

Next, we find that for both U and rj, it is useful to replace the given variable by 
that which involves the subtraction of the indicated log term in the asymptotic VTD 
expansions Eq. (3.8)-(3.13): We set ff := r) - \{1 - k) 2 logt and set U := U - \{1 - 
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k(x))logt; compare this to our approach in Section 2.4.4. Adding a few other minor 
modifications, we are led to define the following set of first-order variables: 





= u, 


u 2 = DU, 


u 3 = td x U, 


(3.17) 


U4 


= A, 


u 5 = DA, 


u 6 = td x A, 


(3.18) 


u 7 


= V, 


u 8 = a, 


u 9 = £. 


(3.19) 


at 


this stage, k(x) 


is an arbitrary function 


(introduced in Eqs. 


(3.8)- 



(3.13)), with no restrictions. In terms of the new set of the variables, the main evolution 
system Eqs. (3.1), (3.2), (3.4) and (3.6) can be written in symmetric hyperbolic form as 
follows: 

Dui-u 2 =0, (3.20) 

Du 2 - u 8 td x u 3 = -tu 9 (u 3 - -t log tk') + ^e 4u H~ 2k (u 2 5 - u 8 u§) (3.21) 

_ \ e lu H l/2{l-k?-2 u%K 2 {l _ k + 2U2) 
--t 2 log tk" U 8 , 

u 8 Du 3 - u 8 td x u 2 - u 8 u 3 = 0, (3.22) 

£>u 4 -u 5 =0, (3.23) 

Du§ — 2ku§ — u 8 td x UQ = — 4u^u 2 + -tugu^ + 2u 8 u§{2u 3 — tlogtk') (3.24) 

_ ^2^1/2(1-^-2^^ 

u s Duq - u$td x u 5 - u$uq = 0, (3.25) 

Du 7 = (1 - k)u 2 + u 2 2 + -u s (2u 3 - t log tk') 2 (3.26) 



or equivalently as 
where 



+ l t ~ 2k e^ {u 2 + u 8 u 2 ) + ^H^-^usK 2 , 
Dug = _ e 2«r t l/2(l-^-2^2 ; (3 _ 27) 

Du9 = _ e 2u 7t l/2(l-kf-2 U8K 2 (3 _ 28) 

f (l - k + 2u 2 ){2u 3 -t\ogtk')u 8 x _ 2k 4 \ 

y Vt u 5 u 6 u 8 e 1 +u 9 J, 

S 1 Du + S 2 td x u + Nu = f[u], (3.29) 
=Diag(l,l,w 8 , 1,1, «s, 1,1,1), (3.30) 



58 



S 2 (u) 



N(u) 



/n 


n 

U 


n 
U 


U 


U 


U 


n 
U 


n 

u 


\ 








— Ug 























-u s 

































































-Us 























-u 8 




































































\o 























o) 


/n 


1 

— 1 


n 
u 


n 
u 


n 
u 


n 
u 


n 
u 


n 
u 


\ 


u 


n 
U 


n 
u 


n 
u 


n 
u 


n 
u 


n 
u 


n 
u 


U 








-u 8 
































-1 


























-2k 





























-Us 

































































Vo 























0/ 



(3.31) 



(3.32) 



Note that we have multiplied the third and sixth equations by Ug. The source-term 
vector / is given by the right-hand sides of the evolution system Eqs. (3.20)— (3.28). The 
reason for keeping this particular form of the matrix N(u) (and not absorbing some of 
its entries into the source-term) becomes clear shortly. 



AVTD solutions of the main evolution system. We now show as an application 
of Theorem 2.4, and as a step towards proving Theorem 3.1, that there exist unique 
solutions to the singular initial value problem for the main evolution system (3.29)- 
(3.32), with AVTD leading-order term 

^0 =("1,0, "2,0, "3,0, "4,0, "5,0, "6,0, "7,0, "8,0, "9,o) 

( I 2k 2k \ (3.33) 

Although not needed for our present argument, we note (by inspecting Eq. (2.46)) that 
this choice of uq is an ODE-leading-order term; cf. Section 2.4.4. 

To check that we have a quasilinear symmetric hyperbolic system, we need to specify 
an exponent vector along with the PDE system and a leading order term. Looking ahead 
to the conditions of block diagonality, we choose 

fi = (fix, Mi, m, (12, ^2,^2, ^3,^4,1^4), (3.34) 

and expect to construct remainders in spaces Xs^ iq with fi given by 

Mi>M3,A t 4>0, [i 2 > 2k. 
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We then find, after replacing u§ by a* + Wg, that so long as we choose a* > 0, and so 
long as we require that all of the asymptotic data functions be contained in some H q (T 1 ) 
(which we fix below), we indeed have a quasilinear symmetric hyperbolic system, which 
in addition does satisfy the block diagonality condition. 

Before continuing the argument that the hypothesis of Theorem 2.4 is satisfied, we 
state our result. 

Proposition 3.2. For any twist constant K G M, for any Sobolev differentiability index 
q > 3, and for any choice of the asymptotic data functions such that A* is an arbitrary 
constant, a*(x) > ; k,U**,a* G H^+^T 1 ), A„ G H q+l (T l ) and 77* G H q {T l ), and k 
satisfies (at each x G T 1 ) either 

(i) k(x) > 1 + y/G (for arbitrary A** the half-polarized case), 

(ii) k(x) > 1 + \/6 or k(x) < 1 — \/Q (for A** = the polarized case), 

there exists a 5\ G (0,<5] ; and a unique solution of the first order main evolution system 
Eqs. (3.29) -(3.32) with leading-order term uq and remainder w G „,q (and Dw G 
Xsx,uq— 1) so long as the exponent vector [i given by Eq. (3.34) satisfies the following 
inequalities for all x G T 1 ; 

1 < m(x) < min {2, (k(x) - 3)(fc(s) + l)/2}, 

- (2k(x) + Vl + 4/c(x) 2 ) < n 2 {x) < 1 + 2k(x), 

< fj, 3 {x) < Hi(x), 

< fx 4 (x) < ~{k(x) - 3)(k(x) + 1). 

Observe here that the inequality just stated for 112 is not required to hold in the 
case of a polarized solution, since in that case A is not a dynamical variable, and this 
condition is vacuous. Although here and below we list results for the polarized and 
half-polarized cases together for compactness, the reader focusing on the polarized case 
may ignore all references to \ii and to W4, w§ and wq. 

As noted above, this proposition is an application of Theorem 2.4 to Eqs. (3.29)- 
(3.32). In the next lemma we verify that under the assumptions of Proposition 3.2 the 
Condition (i) of Theorem 2.4 is satisfied. The first condition follows directly from the 
definition of the energy dissipation matrix Mq. 

Lemma 3.3. The energy dissipation matrix Mq defined in Eq. (2.6) corresponding to 
Eqs. (3.29) -(3.32) , to the leading-order term uq given by Eq. (3.33) and to the exponent 
vector 11 of the form Eq. (3.34) is positive definite at every x, provided that 

o*(x) > 0, H\{x) > 1, 

fj, 2 (x) > max{l, k(x) + + Ak(x) 2 }, /i 3 (x), /i 4 (x) > 0, 

hold for all x G T . 
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The next lemma establishes Conditions (ii) and (hi) of Theorem 2.4. 

Lemma 3.4. The operator 3~(uq) corresponding to Eqs. (3.29) -(3.32), to the leading- 
order term uq given by Eq. (3.33) , and to the exponent vector fi of the form Eq. (3.34) 
satisfies Condition (ii) and (Hi) of Theorem 2.4 for some exponent vector v > fx, for 
some sufficiently small 5 > 0, and for a choice of the differentiability index q > 3, so 
long as a* and 77* are functions in .ff^T 1 ), A** is contained in H q+l (T\), k and [/** are 
elements of H q+2 (T 1 ) , and if at each point x 6 T , the following inequalities hold for 
and k: 

max{0, 1 - (k(x) - 3)(k(x) + l)/2} < m(x) < min{2, (k(x) - 3)(Jfe(x) + l)/2}, 
2k(x) < H2(x) < min{l + 2k(x), ii\{x) + 2k(x)}, 

< /j, 3 (x) < m(x), 

< n4x) < mm{(k(x) - 3)(fc(x) + l)/2,m(x) - 1 + (fc(s) - 3)(k(x) + l)/2}, 
and 

3 < k(x) in the half-polarized case, 
3 < k(x) or k(x) < — 1 in the polarized case. 

In both the polarized and the half-polarized cases, it follows from the two inequalities 
stated above for ji\ that k(x) must either satisfy k(x) > 1 + y/b or k(x) < 1 — y/b. 

Proof. If the operator 3~(wo), defined in Eq. (2.7), is written out explicitly, it consists 
of products of asymptotic data functions, and components of the unknown function w 
(or products involving exponential functions of these) . All of the multiplicands in these 
products are, by hypothesis, contained in designated function spaces (of the form Xs^ q ). 
Thus, to check Condition (ii) of Theorem 2.4, we primarily need to know the multipli- 
cation algebra of spaces such as Xs^q- The result we need is provided by Lemma B.l in 
the appendix. To check Condition (iii), we need results concerning Lipschitz properties 
of products and exponential functions of elements of the spaces Xs^q. Lemma B.3 and 
Lemma B.5 provide these needed results. □ 

Proof of Proposition 3.2. If we wish to use Theorem 2.4 to show that the system dis- 
cussed in Proposition 3.2 admits solutions with the stated properties, it is sufficient that 
i) the asymptotic data functions, which appear in the leading-order matrices 51,0, £2,0 
and A^"o, (i.e. the functions a* and k), be contained in i/' 3+2 (T 1 ) (with q > 3); and ii) 
we choose the function k(x) so that the hypotheses of both of the above lemmas are 
satisfied. We readily check that exponent functions fi±, {12, fJ.3 and which satisfy the 
combined inequalities, can be found if and only if k(x) > 1 + v6 in the half-polarized 
case, and either k(x) > 1 + y/6 or k{x) < 1 — y/E in the polarized case. □ 

The full set of Einstein's vacuum field equations. Thus far, we have constructed 
solutions u of the main evolution equations for the T 2 ~symmetric system with the 



61 



leading-order behavior Eq. (3.33), according to Proposition 3.2. Given such a solu- 
tion u, we may ask under what conditions is this a solution of the full set of Ein- 
stein's vacuum field equations, Eqs. (3.1)-(3.7), with U = u\ + i(l — fc)logt, yl = M4, 
r] = U7 + \{\ — k) 2 log i, and a = ug. 

Proposition 3.5. For any solution of Proposition 3.2 with asymptotic data satisfy- 
ing all the conditions in Theorem 3.1, the full set of Einstein's vacuum field equations 
Eqs. (3.1) - Eq. (3.6) are satisfied, and Eqs. (3.7) can be solved for G\ and G 2 as stated 
in Theorem 3.1. 

Proof. Since the equations for G\ and G2, Eq. (3.7), are semi-decoupled from the rest, 
we ignore them (as well as G\ and G2) to start, and focus on the subsystem Eqs. (3.1)- 
(3.6). To monitor the extent to which this subsystem is satisfied by fields which satisfy 
the main evolution equations, it is useful to define the following set. Based on Eq. (3.3), 
we define 

H := - Vtt +ar lxx + — + — - — + — U t +aU x +^(A t -aA x ) ^-K (3.35) 

and, based on the constraint Eq. (3.5), 

(AU a 

C\ := -?] x + 2tU t U x + —A t A x - (3.36) 

It la 

Based on the constraints which stem from the definition of the new variables which allow 
us to rewrite the original system in first order, we define 

C 2 :=a x - e, (3.37) 

and 

C 3 :=u 2 /t - d t ux, C 4 :=u 3 /t - d x u\, (3.38) 

C 5 :=U5/t - 5 t n 4 , Cq :=u e /t - d x u±. (3.39) 

The constraint- violation terms C3, C4, C5 and Cq are the easiest to handle. Arguing 
as we do in the discussion of A u in Section 2.4.4, we find that the evolution equation for 
each of these, induced by the main evolution system, takes the form -DC3 = 0, DC4 = 0, 
etc. Then, since the form of the leading order term for the main system implies that 
each of these terms must asymptotically vanish, it follows that each must vanish for all 
time. Observe that this determination that 63,0,4, C5, and C% all vanish allows us to 
freely substitute in the consequences of their vanishing in the analysis of H, C\, and 
Such substitution is very useful. 

We now focus on H, C\ and C2, noting that a solution u of Eqs. (3.20) - (3.28) is a 
solution of Eqs. (3.1)-(3.6) (with the above replacements) if and only if H,C\, and C 2 
vanish identically. Presuming that u is a solution of Eqs. (3.20) - (3.28), we calculate 

H = - U8 C M -in 8 , x C 1 + ^(e 4ui ^ 1 " 2 W6 + ^-l-2n 2 )(logtA;'-2^ 1 n3))c' 2 , (3.40) 
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which tells us that so long as we can show that C\ and C2 vanish, it follows that 
H vanishes as well. We may therefore focus on C\ and C2, for which the evolution 
equations take the following form: 

DCl = - }_ K 2 e 2u 7t (k-3)(k+l)/2 ^ 

- [ul - \e iu H~ 2k ul + t log t u 3 k> - )t 2 log 2 t (A;') 2 ) C 2 , 

DC 2 =2e 2u "KH^ k -^ k+1 ^ 2 u 2 d. 

Under our hypothesis (which implies in particular that the coefficients here are contin- 
uous functions of x), this set of evolution equations for C\ and C2 can be treated as an 
essentially independent system of linear homogeneous ODEs at each spatial point. Not- 
ing that (for q > 3) the coefficients on the right side of both equations are well-behaved 
as t — > and converge to zero at every spatial point. Hence this system for the unknowns 
C\ and C2 is of the form Eq. (2.1) with S\ the identity matrix, and S2 and N2 the zero 
matrices. The singular initial value problem of the form C\ = Ci* + w\, C% = C2* + W2 
with w±,w 2 £ ^s,n,g has a unique solution for every prescribed Ci* and C2* and for 
every sufficiently small [i > 0. The definition of the quantities C\ and C2 in terms of 
the variables U, A, 77, ^ and a according to Eqs. (3.36) and (3.37) implies uniqueness 
for all constraint violations which are compatible with solutions of Proposition 3.2. In 
particular, the unique solution corresponding to Ci* = C2* = is C\ = C2 = 0. 

It remains to determine how the condition Ci* = C2* = relates to the choice of 
the asymptotic data functions k, [/**, A*, A**, 7/*, a* and £*. For a solution u as above, 
the functions C\ and C2 defined by Eqs. (3.36) and (3.37) converge uniformly in space 
at t = 0, and we obtain 

C u = ~< + (1 - k)UU - C 2 * = at - 

It follows that C2* = if and only if 

= (3.41) 

and Ci* = if and only if, for an arbitrary constant 770, 

fh(x) = Vo + £ ((1 - k(X))UU(X) - i(loga,)'(X)^ dX. (3.42) 

In particular, the spatial topology implies that we must choose the asymptotic data k, 
[/** and a* such that 

2W ((1 - k(x'))UU(x') - i(loga*) / (x')') da/ = 0. 

We thus conclude that a solution u of Proposition 3.2 is a solution of Eqs. (3.1)- 
(3.6) (with the above replacements) if and only if the asymptotic data functions satisfy 
Eqs. (3.41) and (3.42). 
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It only remains to solve Eqs. (3.7) for G\ and G 2 . The right-hand sides of these 
two equations have the asymptotic behavior (for t — > 0) of a power of t larger than — 1 
uniformly in space. Hence the right-hand sides are integrable in t at t = at every 
spatial point, and the general solution is 

G!(t,x) = G u (x) + ^\ 2v(t '' x) ^Mt F ~L)A(t',x)Kt'- 3 dt', 
Jo 

G 2 (t,x) = G 2 ,{x) - [ e 2 ^rt^/Mj\x~)Kt'-*dt'. 
Jo 

The functions G\ — Gu and G 2 — G 2 * are contained in X^ 1)llhA and Xg litl6tq , respectively, 
for any choice of exponents < ns(x), He(x) < l/2(k(x) — 3)(k(x) + 1), if u is a solution 
of Proposition 3.2. We therefore can take Gu,G 2 * G H q (T ). □ 

3.3.3 Optimal existence and uniqueness result 

The result we prove in Section 3.1 allows for relatively rough asymptotic data, but 
consequently sacrifices some of the expected range (based on numerical and heuristic 
studies [9, 5, 33]) of allowed values for the asymptotic velocity k = k(x). In this section, 
we consider only smooth asymptotic data, and can then prove a result which increases 
the range for k. 

Theorem 3.6 (Optimal result: AVTD (half)-polarized T 2 -symmetric solutions - infi- 
nite differentiability). Suppose one chooses a twist constant K G R, a pair of asymp- 
totic data constants A* and tjq, and a set of asymptotic data functions k,U**,a* (with 
a*(x) > 0), A**,Gu,G 2 * G C 00 ^ 1 ) which satisfy the integrability condition 

J^((l - k(x))UU(x) - i(loga,)'(a:))da: = 0, 

together with, at each point x G T , either 

(i) k(x) > 3 for arbitrary A** (the half-polarized case), or 

(ii) k(x) > 3 or k{x) < — 1 for A** = (the polarized case). 

Then there exists a 5 > 0, and a T 2 -symmetric solution U, A, n, a, G\, G 2 of Einstein's 
vacuum field equations with twist K of the form 

(U,A,r],a,G 1 ,G 2 ) = (U ,A ,r] ,a ,Gi fi ,G 2fi ) + W. 

Here, the leading-order term (Uq, Aq, rjo, ao, Gift, G 2) o) is given by Eqs. (3.8)-(3.13) and 
Eq. (3.14). The remainder W is contained in X$^ tOQ (and DW G Xs^.oo) for some 
exponent vector \x = fx 2 , /m, fJ>5, A*6) with 

Hi,H 2 - 2k, u 3 ,H4, u 5 ,n 6 > 0. (3.43) 

This solution is unique among all solutions with the same leading-order term u$ and 
with remainder W G Xg^oo- 
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Comparing this result with Theorem 3.1, we see the differences in the hypothesized 
regularity of the asymptotic data, and in the allowed range of the asymptotic velocity 
k{x) in the two results. As we find in proving this result, in fact one does not need C°° 
data; data of "sufficiently high differentiability" is enough. One may ask if the reduced 
range for k(x) for rough data is a real effect, or an artifact of our method of proof 
(which remains an open question). Observe one other important difference between 
the two theorems: Theorem 3.6 provides a stronger result regarding the uniqueness of 
solutions to the singular initial value problem for these systems: While in Theorem 3.1, 
there could in principle exist more than one solution for a given set of asymptotic data 
(since /ii has to be larger than one), we find that according to Theorem 3.6, there is 
exactly one solution for the remainder functions in the full space of interest given by 
Eq. (3.43). We note, however, that there may be two solutions for the same asymptotic 
data which differ by a factor which goes to zero faster than every power of t at t = 0. 

The proof of Theorem 3.6 is in many ways similar to that of Theorem 3.1, but with the 
big difference that the latter involves the application of Theorem 2.4, while Theorem 3.6 
is obtained by applying Theorem 2.21. In both cases, these results for Fuchsian singular 
initial value problems (Theorem 2.4 or Theorem 2.21) are applied to the main evolution 
system for T 2 -symmetric solutions. The portion of the proof of both Theorem 3.6 and 
Theorem 3.1 which shows that the existence of a solution for the main evolution system 
implies the existence of a proof to the full system (since one can choose the asymptotic 
data in such a way that the constraints are necessarily satisfied) is essentially the same 
for the two cases. Hence, to prove Theorem 3.6, all we need is the following proposition 
(with the rest of the proof taken care of by the arguments appearing in the proof of 
Proposition 3.5). 

Proposition 3.7. For any twist constant K G R, and for any choice of the asymptotic 
data functions such that A* is an arbitrary constant, a*(x) > 0, k, £7**, a*, t]q G C^iT 1 ), 
and either 

(i) k(x) > 3 (for arbitrary the half-polarized case), or 
(ii) k{x) > 3 or k(x) < — 1 (for A** = the polarized case), 

for all x G T , there exists a 8\ G (0,5], and a unique solution of the main evolution 
system Eqs. (3.29) -(3.32) with leading-order term uq and remainder w G Xsi,u,ao ( an d 
Dw G Xs lt n,oo)> 80 long as the exponent vector fi given by 

(i = (/il,//l,/}l,/i2,M2,^2,M3,M4,M4) 

satisfies the following inequalities for all x G T 1 

< in{x) < min{2, (k(x) - 3)(fc(x) + l)/2}, 
2k(x) < H2(x) < min{l + 2k(x),m(x) + 2k (x)}, 
< /i 3 (x) < m(x), 
< m(x) < (k(x) - 3)(k(x) + l)/2, 

and fi\ strictly smaller than, but sufficiently close to 1 + jjL\ . 
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One of the key differences between the hypothesis here and that of Proposition 3.2 
is the chosen form for the exponent vector \l. The form we use here - in particular, 
the choice of \i\ for the third component - allows for the wider range of k(x). We can 
choose this form here, but not in Proposition 3.2, since here we do only need to satisfy 
the modified block diagonal conditions in Theorem 2.21 (as part of Proposition 2.18). 
Note that we could have carried out a similar modification for the sixth component of 
fi, but it turns out to be unnecessary; the regularity of the spatial derivative of the A 
variable follows automatically. 

Proof. The proof of this result is based on the use of (order n)-leading-order terms, 
which are developed and discussed in Section 2.4. and, in particular, in Theorem 2.21. 
Recalling that the leading order term uq (which takes the form uq with components 
given by Eq. (3.33)) is in fact an ODE-leading-order term, we check the conditions in 
the hypothesis of Theorem 2.21, noting that a portion of these conditions appear in the 
hypothesis of Proposition 2.18. We check that the matrix S^q(uo)Nq(uo) is in Jordan 
normal form, and that our choice of is strictly larger than the negatives of the cor- 
responding diagonal elements of S^qNq. Then, Condition (ii) of Proposition 2.18 (as 
part of Theorem 2.21), can be checked in essentially the same way as is done in the 
proof of Lemma 3.4, but now with the inequalities listed in the hypothesis of Proposi- 
tion 3.7. Condition (iii) of Proposition 2.18 (as part of Theorem 2.21) follows by repeated 
applications of Lemma B.4. □ 



4 Concluding remarks 

Our results here show that there is a large collection of smooth, polarized and half- 
polarized T 2 -symmetric solutions of the Einstein vacuum equations which exhibit AVTD 
behavior in a neighborhood of their singularities. What can we show further? 

Numerical and heuristic studies of T 2 -symmetric solutions [9, 5, 33] strongly indicate 
that AVTD behavior is not found in such spacetimes unless they satisfy a polarization 
condition. These studies do support the conjecture that AVTD behavior occurs generi- 
cally in polarized T 2 -symmetric solutions. While Fuchsian methods of the sort developed 
here are not expected to be effective in determining such genericity, further numerical 
explorations of the polarized T 2 -symmetric solutions could be very useful. Among the 
issues which might be explored numerically is whether the distinction in the results we 
have obtained for solutions of finite differentiability and those which are C°° is significant 
in any sense. 

Observe that one distinguishing feature of our approach here is that an approxima- 
tion scheme is at the core of the method. This scheme can be implemented for numerical 
computations straightforwardly and contains useful built-in convergence and error es- 
timates. In earlier work [12, 14, 15], building on [3], we implemented this scheme in 
the context of semilinear symmetric hyperbolic Fuchsian equations of second-order and 
have obtained very accurate simulations for Gowdy solutions. We expect to get similarly 
good results for the polarized T 2 -symmetric solutions, which we have studied here. Of 
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particular interest would be to explore the issue of whether the "optimal domain" for the 
asymptotic velocity k(x) > 3 (or k(x) < — 1) can only be obtained for smooth solutions 
(as suggested by our discussion in Section 3) and to see what might happen if we try to 
construct such a solution with lower differentiability. 

It is expected, based on numerical simulations [10], that polarized (and half-polarized) 
f/(l)-symmetric solutions exhibit AVTD behavior. Moreover, Fuchsian methods [26] 
confirm this, at least for analytic solutions. The methods we have developed here, gen- 
eralized to PDEs on higher dimensional manifolds (this is done for T n in [2]), should be 
applicable to the polarized and half-polarized ?7(l)-symmetric solutions, showing that 
smooth solutions of this type also exhibit AVTD behavior. 
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Appendices 

A Properties of the spaces X§. 



In this section we list further basic properties of the spaces Xg tfliq which are defined in 
Section 2.2 as the completion of the normed vector spaces (C°°((0, 6] x T 1 ), || • Hs,^), cf. 
Eq. (2.3). Recall that 5 > is a constant, ji is an exponent vector and q is a non-negative 
integer. We now also define the spaces Xg^ q as the set of maps / : (0,5] — > i/ 9 (T 1 ) 
with the property that TZ\p]f is bounded and continuous; cf., Eq. (2.2). If we endow 
Xs^q with the norm || • \\s,n,q, then X$^ q are Banach spaces. Note that if / € A^+e,? 
for some e > 0, then TZ[/j]f : (0, 5] — > H q (T 1 ) is uniformly continuous. 

By definition, all functions in Xg^q can be approximated by smooth functions. Func- 
tions in X§nq, however, can be approximated by a particularly useful sequence of smooth 
functions as follows. 

Lemma A.l. Let f £ Xg^ <q ; i.e., TZ[fi]f : (0,5] — > H q (T l ) is bounded and continuous. 
Let f be defined as follows 



fit) 



f(t), te(0,S], 



Let cj) : R — > R be smooth with <j)(x) > for all \x\ < 1 and <p(x) = for all \x\ > 1, with 
J^4>(x)dx = 1. Let (pLi) be a sequence of positive numbers with limit 0. For any integers 
i,j, we set 

(KMfiijfrx) := [ (nW)(s,y)-<p(^^) —<p(—) dyds. (A.l) 

Jo JtI a i V a i ) a j V a 3 ) 

Then (lZ[fj]f)ij has the following properties: 

(i) (K[Li]f)ij G C°°((0,<5] x T 1 ) for all integers 

(ii) The function 

f ir .= n[^r\nw) hj (A.2) 

has the property that 

fij e X S)IX) q n Xs, m for all integers i,j. 
In particular, for any given integers i,j, one has 

ll(W/)ij(V)||tf«(ri) < C||/ll^, g , for allte (0,6], 
for a constant C > independent oft (but possibly dependent on i, j). 
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(Hi) (Tl[fj]f)i t j(t,x) — > 1Z[fi]f(t, x) fori,j — > oo at a.e. (t,x) G (0,5] x T . 

(iv) If f is such that TZ[fj]f : (0,5] — > ^(T 1 ) is a uniformly continuous map (e.g., if 
f G Xs,[j,+e,q for some e > 0), then 



Wfij - fh 



for i,j — > oo. 



Proof. Observe that TZ[fi]f is a bounded continuous map (0, oo) — > H q (T l ) since TZ[fi] (t)f(t) 
U[ti](5)f(5) for all t > 5. We obtain / G X^^g and ||/||oo, M , g = The first two 

properties of the lemma can be proven by standard arguments. The third one follows 
from Lebesgue's Differentiation Theorem. We only discuss the fourth property. If we fix 
any t G (0,5], then 1 we calculate 

\\K\pL] (IRM-^M/Jy^) " /(*,*)) WhUti) 

= ||(W/)ij(<^)-(W/)(<,x)||fl2(Ti) 



JT 



Cii 



x — y\ 1 
— ' 

Oii J Ctj 



s - t 



a; 



dy ds 



HUT 1 ) 



as a consequence of the condition that f^(j)(x)dx = 1. Now we write 
and therefore obtain the estimate 



(|W- 1 (W/)i J (t,x) - f(t, x )) 



< 



JT 1 

oo 







(W/)(*,y)-(W/)(s,x) )-<j> 



(^M/)( S ,x)-(^[/i]/)(t,x) -0 

/ OLi 



x — y^ 




( s - t 


ai j 




V «i 


x-y\ 




( s - 1 


OLi J 


Ctj 





dy ds 
dy ds 



(A.3) 

HUT 1 ) 
HUT 1 ) 



Writing the first term on the right hand side of Eq. (A.3) as /, we estimate 



I < 



(^M/)(*,y)-(^M/)(«,x)) -<t> 



x-y 



a; 



dy 



1 



hut 1 ) a i 

Now, it is a standard result for mollifiers that for every s G (0, oo) 



s - t 



ds. 



(TZW)(s,y) - (TZW)(s,x)) -cp 

ai 



x-y 
ai 



dy 



HUT 1 ) 



where lim^oo gi(s) = at every s, and for every integer i, the function g is continuous. 
Since TZ\p]f is uniformly continuous, this function gi extends to the interval [0, oo) with 



We write HUT 1 ) to stress that this is the Sobolev space with respect to integration over the x- 
variable. 
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the same properties. Since lZ[fi](t)f (t) = 1Z[fj](5)f (5) for all t > 6, it follows that there 
is a sequence (jji) with limit 0, such that g%(s) < % for all s G [0, oo). Consequently, / 
can be estimated by a sequence (aj, which (i) is independent of j, (ii) is independent of 
t, and (hi) goes to zero in the limit i — > oo. 

We now discuss the second term of the right hand side of Eq. (A. 3), which we label 
as II. The integral over y is trivial, so consequently 



II < 



(K[»}f)(s,x)-(K[»]f)(t,x] 



1 

HUT 1 ) OLj 



ds. 



The term involving the .ff 9 -norm is a uniformly continuous function in s and t. Hence, 
from Lebesgue's Differentiation Theorem and the definition of /, it follows that the s- 
integral converges to for j — > oo, independently of t G (0, 5] and i. This completes the 
proof of the fourth property. □ 

We can now use Lemma A.l to relate the spaces A^g and Xg^q. 

Lemma A. 2. Fix a constant 5 > 0, an exponent vector fx, and a non-negative integer 
q; then for all e > 0, 

Xs,ii+e,q C Xs^^q C X$^^ q . 

Proof. The inclusion Xs^ >q C X§ „q follows easily from the fact that each element in 
Xs,fi,q is the limit of a Cauchy sequence in (C°°((0, 5] x T 1 ), || • H^g), whose elements are 
in particular bounded continuous maps (0,5] —> H q (T l ), and the convergence is uniform 
in time. 

To check the inclusion X$^ +e ^ q C A^ q, let a function / be given in Xs :fJ ,+ e ,q- Hence 
/ satisfies the condition of the previous lemma, in particular that of Condition (iv). It 
follows that / G A^g. □ 

We also wish to comment on time derivatives of functions in A^^ and A^ A . Let / £ 
Xs,^i,q- We say that / is differentiable in time t if the (bounded continuous) map H[fj]f : 
(0,(5] H^T 1 ) is differentiable in the sense of a map between Banach spaces (Prechet 
derivatives). Its time derivative (multiplied by t) D(TZ[fi]f) can then be considered to 
be a map (0,5} -> H^T 1 ), and we set Df := n\}i}- x (D(K\}i}f) - DK\p\f). If this map 
is continuous, then we call / continuously differentiable in t. If this is the case for / and 
if in addition lZ\p]Df is bounded, then we have Df G X$ 4hq . 

Now, let / G Xgnq be continuously differentiable. Then Df is the distributional 
time derivative of / in the following sense. Let (j) be any test function with the 
properties as in Section 2.3.2. Choose e > 0. Then we clearly have that 

s 

)dt = -e (K\p]f, 
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Hence, the boundary term vanishes in the limit e — > 0. The following integrals are 
meaningful for e = 0, and hence we obtain 

! (n\p]Df,cf>) I?{Tl) dt 

,5 ( A - 4 ) 

= -J (pl\p\f, D4>) L2{T1) + (KW + DK[fj]f, 0) i2(T i } ) dt. 
The reader should compare this with the expressions for weak solutions in Section 2.3.2. 

B On products of functions 

We readily check the following results which are useful in dealing with products of 
functions and their relationship to the function spaces X$nq- 

Lemma B.l. Let f G X$^ iq and g 6 Xg„ 2iq be two functions (0, 5] x T 1 — > R, for some 
constant 5 > 0, some smooth exponents [i\ and and an integer q > 1. Then f ■ g is 
in Xs ilXl+IJi2iq and, for some constant C > 0, 

11/ • 5 , II<5, M1 +M2,g < C|l/ll<5, Mi>9 ' Il5l|5,ft2,<?- 
Observe that the condition q > 1 (for one spatial dimension) is essential here. 

Proof. An essential part of the proof of this lemma is the general estimate 

\\f-g\\m<C(\\f\\ m \\g\\ L ™ + \\g\\ H 4f\\ L ~), 

for arbitrary functions / and g in H q n L°°; see Proposition 3.7 in Chapter 13 of [40]. 
The Sobolev inequalities for q > 1 in one spatial dimension then imply 

\\f-g\\m<C(\\f\\ Hq \\g\\ Hq ). 

Working with this inequality, we see that the lemma follows easily if we choose a se- 
quence (fi) which converges to / in the function space Xg^q, and a sequence (ft) which 
converges to g in Xs^ 2 , q , and then write 

fi- 9i- f ■ 9 = fi - {g% - g) + g ■ (fi - /)• 

□ 

Another important result is the following. 

Lemma B.2. Let w be a d-vector-valued function in X$„ q for some exponent d-vector 
H, a constant 5 > 0, and an integer q > 1. Let S be a d x d-matrix-valued function so 
that 1Z[/j] ■ S ■ H[—fi] is an element of X§£ q for an exponent d x d-matrix £ of the form 
£ij = d where £ is an exponent d-vector. Then, the d-vector-valued function S w is in 
x S,(+fi, q and 

\\S W\\sx+H,q ^ CII^M - S - ^[-lA\\s^,q\\w\\s,fJ.,qi 

for some constant C > 0. 
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This lemma is proved essentially in the same way as Lemma B.l. 

Lemma B.3. Suppose that 5 > 0, s > and r > are constants, n, d and q integers 
with d > 1 and q > 1, /i an exponent d-vector, and v\ and U2 exponent scalars. Let 
functions g± , <?2 : U — > M be given where U is an open subset of~M, d . Suppose that 
gi maps all functions w : (0,6] x T 1 — > M. d in B$^^ q>s to elements g\(w) in Bs yUl ,q,r- 
Moreover suppose that there is a constant C\ > with 

\\gi[wi] - ffl[^2]||<5^i,g < CtWw! - W 2 \\ S>m , 

for all w\,W2 : (0, 5] x T 1 — > M. d in -B<5 iAti(?iS . Let us also assume that g2 maps all functions 
w : (0, 5] x T 1 — > M, d in B$^ A ^ S to elements g2(w) in Bs^ 2 ,q,r o- n d that there is a constant 
C 2 > with 

\\g2[wi) - g 2 [w 2 ] 1 1 8,u 2 , q < C 2 \\wi - w 2 \\ s>m , 

for all w\,W2 ■ (0, S] x T 1 — > R d in Bs :ll: q : s- Now, consider h := g\ ■ g 2 , h(w). Then, 
there exists a p > (which is smaller the smaller r is) so that h maps all functions 
w : (0,(5] x T 1 — ?• M. d in Bs^^^ to elements h(w) in B$ Vl+l , 2 ^ p . Moreover, there exists 
a constant C > with 

\\h[wi] - h[w 2 ] \\s, Vl +u 2 ,q < C\\lVl - W 2 \\s tl J,,q, 

for all w\,W2 : (0, 5] x T 1 — > R d in B$^ q ^ s . 

Proof. If w G -B,5 iA t >(?i s, then g±(w) G B^ Vl ^ q ^ and g2{w) G -B<5,i/ 2 ,<?,r- Lemma B.l implies 
that h(w) = gi(w)g 2 (w) G Xg^+^g and 

11^(^)115,^+^2,9 < C \\9l(w)\\s,ux,q\\92(w)\\ S ,u2,q < Cr 2 , 

where C > is the constant in Lemma B.l. This allows us to set p = Cr 2 and hence 
establishes that h(w) G B§ vx + V2q p . Regarding the Lipschitz estimate, we find 

\\t-^ + ^(h[ Wl ](t)-h[w 2 ](t))\\m <C\\t-^(g 1 [w 1 ](t)-g 1 [w 2 ]m\H4t~ U2 92[wi}(t)\\m 

+ C\\t-^ gi [w2}(t)\\ H *\\t- U2 (g2lwi}(t) - g2[w2}(t))\\m- 

(B.l) 

Then we can use the individual Lipschitz estimates for g\ and g2 in order to establish 
this result. 

□ 

While Lemma B.3 is adequate for the proof of Theorem 2.4, to prove Theorem 2.21 
we require a stronger result, which we present here. 

Lemma B.4. Suppose that q > 1. Let g\ and g2 be functions satisfying all the conditions 
of Lemma B.3 with exponents v\,v 2 for all x G T . Suppose that, in addition, one has 
the following: For all w G Bg^ q s / 2 with uj G Bg^ q>s / 2 for some exponent vector p which 
satisfies p> p, there exist scalar exponents 71,72, independent of fx, such that 

gi{w + to) - 01 (to) G X 5)P+7l)? , g 2 {w + u) - g 2 (w) G X SiJ1+1 ^ q , 
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and 

\\gx[w + u] - fl'iMlla^+T!,, < Ci\\u\\ S)%q , 
\\g 2 [w + u] -52MIU,iK+7 2 ,g < C 2 \\uj\\s^, q , 

for constants C\,C 2 > 0. 

Then the function h = g\ ■ g 2 has the following property. We can choose a scalar 
exponent 7 smaller or equal than min{z^i + 72, v 2 + 71} (independently of fi), such that 
for all w G B$^ qs/2 and u G B 8 ^ q!S / 2 , one has 

h(w + u) - h(w) G X 5i £ +7i9 , 

and 

\\h[w + u\ -h[w]\\ sfi+ ^ q < C\\u\\ S fc q 

for a constant C > 0. 

This follows from a more detailed analysis of Eq. (B.l). 

To handle the exponential function, we rely on the following result. 

Lemma B.5. Pick constants 5 > 0, s > 0, an integer q > 1, and an exponent (i > 0. 
Let := expollj, where II, : ~R d — > M is the projection to the ith component of d- 
vectors. Then, for every function w : (0, 5] x T 1 — > M in B^^^ )S , there exists an r > 0, 
so that the composed function o w : (0, 5] x T 1 — > K zs in i?<5,o,g,r- Moreover, for all 
wi,w 2 G B$^^ s , there exists a constant C > 0, so that 

\\9®(wi) ~ 9 {l) (W2)\\8,0,q < C\\WI - W 2 \\ 5 ^ q . 

In addition, for every scalar exponent fi > /i and every w G B 8liqs / 2 and u G Bg^g^, 
it follows that g^ (w + uj) — (w) is in Xg^q and the estimate 

\\g {i) (w + uj) - g {t) (w)\\ SAq < C\\uj\\s Aq , 

holds. 

Proof. This follows from Proposition 3.9 in Chapter 13 of [40] applied to g^(w) — 1, 
together with the Taylor theorem for the exponential function. □ 



C Duality and convergence results 
Sobolev spaces and duality 

Following [16, Chapter VI] or [36], one defines the Sobolev space H s (M. n ) for any s£R 
as the set of temperate distributions u such that u(l + |£| 2 ) 5 / 2 G L 2 (M n ), where u := Tu 
is the Fourier transform (in the sense of temperate distributions) of u. The norm defined 
by 

\\u\\ s := + |£| 2 r /2 |lL ? (R") 
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turns this space into a Banach space. If s = q for any non-negative integer q, then 
H s (R n ) is equivalent to the standard (p = 2) Sobolev space H q (R n ). For general s£l, 
the space H s (R n ) is in fact a Hilbert space for the scalar product 

(u,v) s := [ «(0(1 + |£| 2 r /2 £(6(l + \t\ 2 ) s/2 dt. 

Let u G // _s (E n ) and u G .£P(]R n ) for any s£i Then the dual pairing between 
# s (IR n ) and ii^ s (M n ), 

(«,*):= / ( C1 ) 

JR™ 

is well-defined, as a consequence of the inequality 



(u,v) | < 



^)(1 + IC| 2 )- S/2 ^)(1 + IC| 2 ) S/2 ^ 



< H_ s ||u|| s . (C.2) 



By means of this pairing, we can identify H~ s (R n ) with H s (R n )* (the dual space) 
as follows. For every u G H~ s (R n ), the map (u, •) : H s (R n ) — > R is a bounded linear 
functional, i.e., an element of i/ s (]R n )*. Conversely, according to the Riesz representation 
theorem, there exists a unique element G H s (R n ) for each element (f> G fP(R n )* such 
that 

4>(v) = {w 4> ,v) s 
for all v G -fP(]R n ). The last expression can be written as 

where := w$(tj)(l + |£| 2 ) s is the Fourier transform of := J r_1 ( , u; <) i(^)(l + |£| 2 ) s ). We 
have V4> G H- s (R n ), since %(1 + |£| 2 )" s / 2 = + |£| 2 ) s / 2 G L 2 (U). By means of the 

pairing above, we have thus constructed a unique element G H~ s (R n ) corresponding 
to each cfi G H s (R n )*. In this sense, we can therefore identify H~ s (R n ) with H s (R n )* 
for every s£R. 

The following result concerns the relationship between Sobolev spaces of different 
indices. 

Proposition C.l. For every s£l and a > 0, the space H s+a (R n ) is a dense subset of 
H s {R n ). 

Proof. We first show that H s+a (R n ) is indeed a subset of H s (R n ) for a > 0. Suppose 
that tt G -fP +<J (R n ). Calculating the || • || s norm of u, we obtain 

imi 2 = / is(e)i 2 (i + ici 2 m< / is(oi 2 (i + ici 2 r +,T ^ = iiuii 2 +CT <oo, 

Jr™ Jr 71 

from which it follows that u G _fP(M n ). To check that H s+(T {R n ) is a dense subset, it 
is sufficient to note (see, e.g., [16]) that C^°(R n ) (the space of smooth functions with 
compact support) is dense in both H s (R n ) and H s+a (R n ). □ 
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Convergence results in Sobolev spaces 

One can use this dense inclusion property (Proposition C.2) together with the duality 
properties discussed above to derive certain convergence and closedness-type results for 
sequences in Sobolev spaces. We first discuss a result of this sort for Sobolev spaces on 
M n , and then do the same for Sobolev spaces on T 1 . 

Proposition C.2. Choose s,sq G R so that < so < s. Let (w m ) be a bounded 
sequence in H s (M. n ) in the sense that there exists a constant C > so that \\w m \\ s < C, 
for all integer m. Moreover, suppose that (w m ) converges to some w G H S0 (M. n ); i.e., 
\\w-m — w \\s ~ > 0- Then, w is contained in H s (M. n ). 

Proof. The boundedness of the sequence implies the existence of a subsequence of (w m ) 
(which for simplicity we identify with (w m )) which converges weakly. Hence, as a con- 
sequence of the Riesz Representation Theorem and the above dual pairing in Eq. (C.l), 
there exists an element w G H s (M. n ), so that, for every Y € H~ s (M. n ), 

(Y,w-w m )^0 (C.3) 

We wish to show that w = w and hence that w € H s (M. n ). To do this, we consider an 
arbitrary X G H~ s °(M. n ) and the dual pairing 

\(X,w- w)\ <\(X,w- w m )\ + \(X,w- w m )\ , 

where w — w is considered as an element of H~ s °(K n ), and where we have used the 
triangle inequality. Since X £ H~ s °(M. n ) C H~ s (M. n ) according to Proposition C.l, we 
can consider the first term on the right hand side as a pairing between H s (M n ) and 
H~ s (M. n ), and hence Eq. (C.3) implies that this term can be made arbitrarily small 
by choosing m sufficiently large. The second term is considered as a pairing between 
H S0 (R n ) and H~ s °(R n ) so that Eq. (C.2) yields 

\(X,w-w m )\ < ||X|| 

— so\\ w w rn\\s - 

Also this term can be made arbitrarily small by choosing m sufficiently large. Hence, we 
have found that (X,w — w) = for all X G H~ s °(M. n ). Now, the Riesz representation 
theorem implies that for every X G -ff _s °(IR n ) there exists precisely one X G H s °(W l ) 
for which 

0=(X,w-w) = (x,w-w) 

\ I H 3 (R n ) 

In particular, we may choose X = w — w, which implies that w — w = 0. □ 

Corollary C.3. Choose non-negative integers q and qo so that qo < q. Let (w m ) be 
a bounded sequence in H q (T l ), in the sense that there exists a constant C > so that 
H^mll/f^T 1 ) — C, for all integers m. Moreover, suppose that (w m ) converges to some 
w G {{^(T 1 ); i.e., \\w m — ■u;||^9 (t 1 ) — * 0- Then, w is contained in i/ 9 (T 1 ). 



77 



Proof. We formulate the proof so that it can be easily generalized to general smooth 
orientable, connected compact Riemannian manifolds M in any dimension n. For this 
paper, the relevant special case is M = T . Let ((Ui,fa)) be a collection of coordinate 
charts, i.e., open subsets Ui C M and homeomorphisms fa : Vi — > U{ where Vi := fa[ (Ui) 
are open subset of M n , which cover M, i.e., M = U i L^. It follows from compactness 
that we can assume that there are N such coordinate charts. Let (r^) be a subordinate 
partition of unity. Then we find that (w m ) is a bounded sequence in H q (M) if and only 
if for all i = 1, . . . , N, we have that (w m ofa) is a bounded sequence in H q (Vi). Moreover, 
\\w m — w\\jjqo (t 1 ) ~~ ^ f° r some w G H qo (M) if and only if for all i = 1, . . . , N, we have 
that \\w m o 4>i — w o 4>i\\mo(Vi) ~^ (since w o fa £ H qa (Vi)). Now, the Stein Extension 
Theorem (Theorem 5.24 in [1]) implies the existence of total extension operators 
Ei (Definition 5.17 in [1]), which are linear maps Ei from functions defined on Vi to 
functions defined on M. n with the following property: If / € H r (Vi) for any non-negative 
integer r, then 

1. (Eif)\y. = f almost everywhere, 

2. Eif is in H r (M n ), and there exists a constant C > 0, so that 

\\Eif\\H r (RK) < C\\f\\ H r(y.y 

Hence, we find that (w m ) is a bounded sequence in H q (M) if and only if for all i = 
1,...,N, we have that (Ei(w m o fa)) is a bounded sequence in H q (M. n ). Moreover, 
\\w m — if 1 1 ao (t 1 ) — > f° r some w G H qo (M) if and only if for all i = 1, . . . , N, we have 
that \\Ei(w m o fa) — Ei(wo fa) \\ H q / R „) — )■ (since Ei(wofa) £ H qo (W 1 )). It follows from 
Proposition C.2, that ^(io o f ) G i7 g (M n ). Hence, w o fa £ H q (Vi). Since this is true 
for alH = 1, . . . , N, it follows that «; G H q (M). □ 
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